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ABSTRACT 


A  new  method  for  absolute  calibration  of  electroacoustic  trans¬ 
ducers  is  described  which  permits  low  frequency  reciprocity  calibration 
in  a  test  apparatus  of  small  dimensions.  The  method  takes  advantage  of 
the  reduced  wave  speed  in  a  cylindrical  water  column  bounded  by  a 
compliant  PVC  tube  to  reduce  the  length  of  the  standing  wave  to 
dimensions  which  can  be  easily  handled  in  the  laboratory. 

-  Results  of  a  limited  series  of  experiments  show  a  mean  reproduci¬ 
bility  of  ±0.8  decibels  in  a  one-meter  long  tube  at  frequencies  between 
750  and  1100  Hz  and  sound  speeds  from  330  to  365  meters/second.  The 
hydrophone  voltage  sensitivities  obtained  by  the  method,  however, 
systematically  differ  by  1.5  decibels  from  those  measured  by  comparison 
with  a  standard  hydrophone;  additional  research  will  be  required  to 
resolve  this  discrepancy. 
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I.  WHY  COMPLIANT  TUBE  RECIPROCITY  WORKS 


A.  ACOUSTIC  CALIBRATION  AND  RECIPROCITY 

There  are  two  general  classes  of  underwater  acoustic  transducer 
calibration:  absolute,  or  primary,  calibration,  and  comparison,  or 

secondary,  calibration.  The  latter  presumes  the  availability  of  a  sound 
field  of  known  amplitude  measured  by  a  calibrated  transducer,  but 
absolute  calibration  can  be  performed  using  only  voltage  measurements 
and  primary  measurements  of  the  test  system  environment  in  a  sound 
field  of  unknown  amplitude.  Once  the  voltage  sensitivity  of  a  standard 
hydrophone  has  been  established  in  a  particular  frequency  range  and 
environment  by  absolute  calibration,  it  may  be  used  as  a  reference  for 
comparison  calibration  of  other  hydrophones  and  of  projectors;  this 
absolute  calibration  is  most  often  done  by  the  reciprocity  method. 

1.  Fundamentals  of  Reciprocity  Calibration 


The  reciprocity  method  for  absolute  calibration  is  based  on  the 
properties  of  linear,  passive  systems  whose  action  can  be  described  by 
the  two-port  equation  (1.1),  which  relates  the  acoustic  pressure  and 
velocity  (p,u)  at  a  transducer's  mechanical  terminal  and  the  voltage  and 
current  (E,l)  at  its  electrical  terminals  by  a  symmetric  or  antisymmetric 
matrix  of  constant  coefficients;  such  a  system  is  called  reciprocal. 


\  _  ’  *  "  / u 

/  L±$  Zr  J  \  | 


(1.1) 


Hunt  [Ref.  1]  derives  this  equation  from  first  principles  and  shows  how 
the  matrix  elements  may  be  evaluated  for  common  transducers.  Foldy 
and  Primakoff  show  [Ref.  2]  that  the  short  circuit  transmitting  current 
response  S  and  the  open  circuit  hydrophone  voltage  response  M  of  a 
reciprocal  transducer  T  are  related  by  a  constant  factor 


-  »  .  -  ,  |  >  .  '  «  .  H  H  ,  «  •  -  •  »  «  «•  \  \  V  ",  %  %  •,  , 


where  J  is  called  the  reciprocity  parameter  and  is  characteristic  of  the 
test  medium  and  its  boundaries.  Bobber  has  demonstrated  [Ref.  3] 
that  J  is  equivalent  to  the  acoustic  transfer  admittance  between  an 
arbitrarily  defined  transmitting  surface  T  generating  a  volume  velocity 
U  and  an  arbitrarily  defined  receiving  surface  H  within  the  test  envi¬ 
ronment  exposed  to  an  average  pressure  P,  and  exists  as  a  function  of 
:he  test  medium  and  its  boundaries  for  any  arbitrary  linear,  passive, 
and  reversible  medium. 

When  two  test  transducers1  are  located  at  points  in  the 
medium  where  the  pressure  amplitudes  are  equal,  the  reciprocity  param¬ 
eter  allows  expression  of  the  signal  voltage  of  a  test  hydrophone  H 
placed  at  one  point  as  a  function  of  the  transmitter  current  response  of 
a  reversible  transducer  T  placed  at  the  other  point  and  driven  to 
create  the  sound  field; 

Eth  =  MhPh  —  MhStIt  (1.3) 

where  I  is  the  drive  current  input  to  T.  If  the  sound  field  generated 
by  T  is  then  replaced  by  a  field  which  is  identical  at  the  locations  of  T 
and  H,  but  is  generated  by  an  separate  projector  P,  the  outputs  of  T 
and  H  in  this  new  field  will  be  related  by  equation  (1.4).  Combining 
this  equation  with  (1.3)  produces  the  calibration  equation  (1.5),  which 
allows  the  computation  of  the  test  hydrophone's  voltage  sensitivity  H 
directly  from  voltage  measurements  made  in  the  three  sou rce-to- receiver 
configurations:  T-to-H,  P-to-H,  and  P-to-T. 

lTo  avoid  confusion,  the  term  test  transducers  will  always  be  used 
when  referring  to  both  transducers  used  in  the  reciprocity,,  cal ibration . 
The  term  reversible  transducer,  or  the  ANSI  designator  "T",  will  be 
used  to  refer  to  the  reciprocal  transducer  used  both  as  a  hydrophone 
and  as  a  projector  during  the  calibration;  test  hydrophone,  or  "H"  will 
designate  the  transducer  used  only  as  a  hydrophone;  and  projector,  or 
’P”  will  identify  a  third  sound  source  used  to  create  the  sound  field 
used  in  the  test. 


M.P  J  S 


Mh  =  J 
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rocity  Test  Geometries 


Most  reciprocity  calibrations  are  made  in  one  of  the  two 
earliest  system  geometries  for  which  the  reciprocity  parameter  was 
computed:  free  field  [Ref.  4],  or  close-coupled  [Ref.  5].  The  tech¬ 
niques  for  carrying  out  these  tests,  and  the  limitations  required  to 
ensure  the  validity  of  the  test  are  described  by  the  ANSI  standard  for 
underwater  calibration  [Ref.  6].  The  free  field  test  geometry  requires 
a  volume  of  water  large  enough  to  permit  spherical  spreading  of  the 
acoustic  waves  undisturbed  by  boundary  reflections.  This  implies  tanks 
of  dimensions  greater  than  ten  wavelengths  (about  15  meters  at  1  kilo¬ 
hertz)  for  underwater  testing.  Coupler  reciprocity  is  limited  by  the 
requirement  to  encase  the  entire  active  surface  of  two  test  transducers 
in  a  rigid-walled  cell  of  dimensions  less  than  1/10  of  a  wavelength; 
useful  underwater  low  frequency  transducers  limited  to  such  dimensions 
can  be  difficult  to  build. 

The  practice  of  reciprocity  calibration  is  not,  however, 
constrained  to  the  two  conventional  geometries.  Rudnick  [Ref.  7] 
discusses  three  compact  geometries  with  known  reciprocity  factors; 
Bobber’s  reference  text  on  measurements  [Ref.  8]  describes  several 
more.  These  unusual  test  environments  extend  the  practice  of  reci¬ 
procity  calibration  beyond  the  limits  of  of  the  two  conventional  geome¬ 
tries.  A  summary  of  these  methods  is  included  in  Table  I;  each  of  the 
reciprocity  parameters  listed  is  the  quotient  of  an  area  characteristic  of 
the  system  geometry,  and  of  the  specific  acoustic  impedance  of  the 
medium. 

The  compliant  tube  test  system  described  in  this  paper  adds 
another  option  to  the  array  already  available  to  the  experimenter .  It 
has  particular  value  to  the  .  underwater  acoustician  because  it 


TABLE  I 

Some  Reciprocity  Calibration  Systems 


Geometry 

Reciprocity  Parameter 

Reference 

Free  field 

2dX 

Po  c 

4 

Close— coupled 

2sAjj£ 

PoC 

5 

Plane— wave 

2A. 

PoC 

8 

Cylindrical— wave 

2LV5X 

P  qO 

8 

Helmholtz  resonator 

27TALf 

Po  c*Q 

7 

Waveguide, 

resonant  plane  wave 

7rALf 

PqcQ 

7 

significantly  lowers  the  frequency  limit  for  underwater  reciprocity  cali¬ 
bration  below  the  practical  limits  of  the  conventional  geometries. 
Unlike  rigid-tube  methods,  it  can  be  accomplished  within  a  convenient, 
inexpensive  test  cell,  and  although  still  limited  by  the  requirement  for 
accurate  measurement  of  the  system  quality  factor,  does  escape  resonant 
reciprocity  calibration's  usual  restriction  to  discrete  frequencies.  The 
"virtual  rigid  walls"  created  by  the  standing  wave  pressure  antinodes, 
and  in  which  the  transducers  are  placed  for  the  test,  can  be  adjusted 
to  any  resonant  length  and  frequency  by  simply  raising  or  lowering  the 
water  level. 

Our  development  of  the  compliant  tube  reciprocity  parameter 
must  begin  from  an  understanding  of  wave  motion  in  the  water  column, 


and  of  the  oscillatory  expansion  of  the  tube's  circumference.  From 
these  equations,  we  will  develop  the  energy  balance  of  the  oscillating 
system;  this  in  turn  will  provide  the  result  for  the  reciprocity 
parameter. 

B.  COMPLIANT  TUBE  WAVE  EQUATION 
1 .  Waves  m  the  Water 

Lighth ill's  development  [Ref.  9]  of  the  wave  equation  for  a 
lossless  fluid  moving  coaxially  in  a  lossless  tube  of  varying  cross 
section  includes  the  case  of  the  compliant  elastic  tube  in  which  our 
experiment  is  performed.  The  momentum  and  mass  conservation  equa¬ 
tions  for  an  inviscid,  ideal  fluid  are 


d(pq) 
a  t 


=  -2e 

dz )  3z 

'  r(pau) 


(1.6) 

(1.7) 


where  p  is  the  acoustic  excess  pressure,  u  the  acoustic  coaxial  particle 
velocity,  and  a  the  local  time-varying  cross-sectional  area  of  the  tube, 
which  varies  slightly  about  its  equilibrium  value  A.  When  fluid  entropy 
is  constant,  both  the  local  density  p  and  the  local  cross  section  a  may 
be  expressed  as  functions  of  the  acoustic  pressure.  (The  latter  relation 
is  just  the  inverse  of  the  function  which  relates  pressure  to  cross 
section  variations.)  Discarding  terms  of  second  order  and  higher 
reduces  these  equations  to 


du 

_^P 

(1.8) 

Poat~ 

dz 

<3(pa)_ 

.du 

-po  A— 

(1.9) 

a  t 

dz 

(1.9) 


Taking  the  spatial  derivative  of  (1.8)  and  the  time  derivative  of  (1.9), 
combining  the  results  with  the  linearized  equation  of  state  p  =  (p  -p9)c2 
and  discarding  terms  above  first  order  gives  a  wave  equation 


d2p  _  -2  d2p 
dz  ~  °  3?" 


where  the  wave  speed  c  is  defined  by 


(1.10) 


-2  __  (j_  d(pa)\  _  /  l  d£ 

P\p0A  d  p  /o  p0  \  p0apj 


i*l)  +  D) 


(1.11) 


(The  subscript  "0"  above  implies  evaluation  of  the  derivative  at  the 
undeformed  wall  position  and  under  adiabatic  conditions.)  The  first 
term  of  equation  (1.11)  is  the  adiabatic  compressibility  or  reciprocal 
bulk  modulus  K;  if  there  were  no  dependence  on  cross  section,  this 
would  reduce  to  the  equation  for  the  sound  speed  in  the  bulk  medium. 
If  the.  second  term,  called  the  distensibility  D,  is  much  larger  than  the 
adiabatic  compressibility,  the  speed  of  sound  waves  in  the  tube  will  be 
determined  by  the  characteristics  of  the  wall  rather  than  by  those  of 
the  bulk  fluid,  and  will  be  much  slower.  This  in  turn  will  reduce  the 
size  of  the  resonant  tube  required  to  operate  at  a  given  frequency. 

The  harmonic  solution  of  equation  (1.10)  with  a  pressure 
release  boundary  at  z=0  and  a  rigid  boundary  at  z=L,  and  for  frequen¬ 
cies  below  the  cutoffs  for  azimuthal  and  radial  modes,  is 


p(z,t)  =  P  sin(kz)sin(cjf) 


(1.12) 


with  a  corresponding  velocity  field 


(M3) 


u(z,t)  = - cos(kz)  cos(«t) 

P  oc 

The  wave  number  k  is  limited  by  the  boundary  conditions  to  the 
discrete  values 


k  = 


27T 

\ 


■n~(n  +  1/2) 
L 


n-0,  I*  2f  3,.... 


(1 ■ 14) 


This  sound  field  has  associated  with  it  both  kinetic  and  potential  energy 
density.  The  kinetic  energy  of  axial  water  motion  in  a  thin  radial  cross 
section  of  length  dz  is 


(1.15) 


The  potential  energy  in  the  same  section  is  found  by  considering  the 
work  done  by  the  fluid  as  the  tube  boundaries  expand: 


Adz 


1  do  1 

AdpJa 


.  .  1  dp] 

-  A  dz - c 

Podpl 


dp 


=  (K  -  D)f 


Adz 


(1.16) 

(1.17) 


where  £  is  a  coordinate  of  the  periodic  motion. 

2.  Motion  of  the  Wall 

The  distensibility  can  be  calculated  from  the  material  proper¬ 
ties  of  the  tube  wall  and  its  geometry.  In  the  case  of  a  cylindrical  tube 
which  responds  elastically  to  small  excess  stresses,  the  circumferential 
tensioning  force  component  due  to  the  pressure  on  a  thin  radial  cross 
section  of  length  dz  is 


0.19) 


<$o  F  _2R 

£  “  A  ~  2hYdz  “  hY  P 

This  tensioning  force  causes  a  strain  s  defined  by  equation  (1.19). 
Where  p  is  small  and  a  varies  only  slightly  from  its  undisturbed  value 
A,  the  distensibility  is  then 


A  dpJQ  hY 


(1.20) 


This  stretching  of  the  wall  along  its  circumference  has  the  characteris¬ 
tics  of  a  damped,  stiffness-controlled  driven  oscillator  (neglecting  long¬ 
itudinal  effects).  The  ratio  of  the  driving  pressure  at  any  cross 
section  and  the  circumferential  strain  velocity  is  the  mechanical  impe¬ 
dance  of  the  cross-sectional  ring  treated  as  a  damped  oscillator  of 
negligible  mass  with  elastic  constant  Y  and  damping  coefficient  R(m) 
[Ref.  10:  pp.  10-12] 


2P  sin(kz)sin(cjt-0) 

s  = - 4 — 1 - 5 - 

cjW 

w  =  yjni  +(Y/u)z  9  =  arctan^^-  j 

The  polyvinyl  chloride  (PVC)  tube  used  in  this  experiment  is  actually 
not  elastic  over  any  finite  range  of  distortion,  so  the  parameter  Y  is 
not  Young's  modulus;  it  is  only  a  linear  approximation  to  the  plastic 
deformation  that  the  PVC  undergoes  in  harmonic  motion.  Within  this 
linearization,  the  potential  energy  stored  in  a  thin  radial  section  of  the 
wall  as  the  tube  is  expanded  by  a  small  amount  from  its  initial  cross 
section  can  be  found  from  the  force-displacement  integral 


(1.21) 

(1.22) 


dE#  =  J  F  ~y~  as  =  Dp  J  Adz 


0.23) 


with  the  force  given  by  equation  0.18) 


C.  COMPLIANT  TUBE  RECIPROCITY  PARAMETER 

The  arguments  used  by  Garrett,  et  a}.,  [Ref.  11]  to  derive  the 
reciprocity  parameter  J  for  a  rigid  tube  will  also  yield  an  expression  for 
the  reciprocity  parameter  of  the  free  surface,  compliant  wall  tube.  As 
with  the  rigid  tube,  the  solution  depends  on  relating  the  energy 
coritent  of  the  system  E  and  the  energy  loss  per  cycle  AE  to  the 
acoustic  transfer  admittance  (reciprocity  parameter)  described  by  equa¬ 
tion  (1.2).  When  the  ratio  of  energy  content  to  energy  loss  rate  is 
known  in  terms  of  J  and  of  primary  measurements,  the  resonance 
quality  factor  equation  (1.24)  (which  can  also  be  found  from  primary 
measurements)  can  be  used  to  evaluate  J. 


f„  -  f,  AE/2tt 


(1.24) 


1 .  Svstem  Enerqv  Content  and  Losses 

The  compliant  tube  system  stores  energy  both  in  the  water 
column  and  in  the  tube  wall.  Kinetic  energy  of  motion  of  the  water 
pressurizes  the  water  locally,  which  in  turn  stresses  and  expands  the 
tube  diameter.  The  energy  balance  exchanges  between  the  three 
degrees  of  freedom  are  governed  by  the  work  equations  for  pressuriza¬ 
tion  (equation  (1.17))  and  for  expansion  (equation  (1.23)),  and  by  the 
kinetic  energy  equation  (1.15).  Summing  these  equations  gives  the 
total  energy  in  a  thin  cross  section 


^ -  -  g ~  s 1  ■.  I',  rg ■-■  U  ^*7  V’  <■ 'J '  '** r* 
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=  Adz  j  +  ^(K  +  0)  |=  Adz  |  ^ 


(K  -  D)|-+  Dp'  | 

(1.25) 

(  2  2p0c  ) 

(1.26) 

where  equation  (1.11)  has  been  used  to  replace  K  and  D  with  p0c2. 
The  sound  field  in  the  absence  of  attenuation  is  described  by  equations 
(1.12)  and  (1.13);  substitution  gives  a  form  which  can  be  integrated 
over  the  length  of  the  water  column 


p2  ) 

+ -  2  sin2(kz)  sin2(ut)£  dz  (1 .27) 

2p0c  ) 


(1.28) 


This  is  identical  to  the  result  for  a  rigid  resonant  tube  [Ref.  11]. 

In  the  steady  state,  the  cyclic  energy  loss  equals  the  work 
done  by  the  driving  transducer.  If  this  transducer  is  located  at  a 
pressure  antinode,  pressure  and  velocity  must  be  in  phase  at  the 
transducer  face  (  since  the  impedance  to  flow  at  a  velocity  node  is  infi¬ 
nite).  The  energy  transferred  to  the  system  is  the  work  done,  which 
can  be  written  as 


n  ds 


=  P  U  sin^wt) 


(1.29) 


where  U  is  the  volume  velocity  amplitude  of  the  transducer  surface. 
Integrating  over  one  cycle  gives 


(1.30) 
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2.  Reciprocity  Parameter 


Our  calculations  above  have  yielded  both  the  system  energy 
content  and  the  system  energy  loss  (actually,  its  throughput  at  steady 
state).  The  energy  equation  for  Q  (1.24)  immediately  produces 

Q  E  _  P  ALf 
2n  AE~U2p0c2 

The  ratio  of  volume  velocity  amplitude  to  pressure  amplitude  that 
appfears  in  this  result  is  the  acoustic  transfer  admittance  between  pres¬ 
sure  maxima  (equation  (1.2)),  and  must  also  be  the  reciprocity  param¬ 
eter  for  hydrophones  positioned  at  separate  maxima  when  the  column  is 
driven  at  resonance.  Equation  (1.32)  follows;  again,  it  is  the  same  as 
the  reciprocity  parameter  for  the  rigid  tube. 


7TALf 

PoCZq 


(1.32) 


The  solution  of  equation  (1.32)  requires  the  experimental  measurement 
of  three  resonance  parameters:  the  resonance  frequency  f;  the  reso¬ 
nance  quality  factor  Q;  and  the  wavelength'  of  the  resonance  standing 
wave,  in  order  to  compute  the  sound  speed  C.  The  accuracy  of  the 
value  computed  for  J  depends  on  our  ability  to  measure  these  parame¬ 
ters  accurately.  The  methods  for  doing  this  will  be  discussed  in  the 
next  chapter. 

In  using  equation  (1.5)  with  the  result  derived  above,  we 
must  ensure  that  the  pressures  and  voltages  implied  are  indeed  the 
maximum  amplitude  values  at  the  antinodes.  The  test  transducers' 
output  will  be  proportional  to  the  average  sound  field  over  their  active 
surfaces,  which  at  its  maximum  will  be  less  than  that  which  would  be 
generated  by  a  uniform  field  at  the  maximum  pressure  amplitude  P  by 
the  factor 


in 


[o: 


Measured  voltage  outputs  must  be  corrected  by  this  factor  before  the 
results  are  used  in  the  calibration  equation  when  transducer  dimensions 
are  an  appreciable  fraction  of  a  wavelength.  For  the  LC-10's  used  in 
our  experiments,  this  factor  is  0.05  dB  at  1kHz;  for  the  LC-32's,  it  is 
0.27  dB. 

3.  Limits  of  Derivation 

The  derivation  of  the  compliant  tube  reciprocity  parameter  has 
been  done  under  four  important  simplifying  assumptions.  To  the  degree 
that  these  approximations  do  not  represent  the  real  system,  the  results 
obtained  by  absolute  calibration  will  be  in  error. 

The  first  of  these  approximations  is  the  use  of  the  solutions 
to  the  lossless  wave  equations  in  computing  the  total  energy  in  equation 
(1.28).  The  PVC  wall  does  attenuate  the  sound  wave,  resulting  in  a 
resonance  sound  field  which  resembles  that  derived  by  Kinsler,  Frey, 
Coppens,  and  Sanders  [Ref.  10: pp. 206-209] .  Figure  1.1  shows  the 
measured  variation  of  the  pressure  amplitude  with  depth  for  the  n=4 
mode.  The  pressure  amplitudes  at  separate  maxima  will  therefore  not 
be  identical,  as  is  required  by  our  development,  which  assumes  an  infi¬ 
nite  standing  wave  ratio. 

A  second  approximation  is  the  assumption  of  linear  behavior 
for  the  polyvinyl  chloride  wall.  In  fact,  PVC  has  no  elastic  region  for 
strain,  as  was  discussed  in  the  derivation  of  the  wall  motion  equation. 
The  nonlinearity  of  the  wall  permits  mode  coupling  into  other  resonance 
modes,  particularly  when  two  resonances  have  spatial  pressure  maxima 
near  the  same  position  and  are  related  subharmonically .  The  result  of 
such  coupling  when  the  coiumn  is  driven  at  one  frequency  is  to  depress 
the  amplitude  of  individual  pressure  maxima  by  an  amount  proportional 
to  the  degree  of  coupling  at  that  position;  again  violating  the  require¬ 
ment  for  equal  pressure  amplitudes  at  the  two  test  positions. 


Third,  the  reciprocity  concept  requires  that  the  receiving 
transducer  approximate  an  infinite  mechanical  impedance  to  the  medium. 
While  piezoceramics  do  meet  this  approximation  very  well,  the  soft 
mounting  material  used  the  the  hydrophone  body  and  gas  bubbles  on 
the  hydrophone  surface  do  not.  Their  presence  could  significantly 
distort  the  sound  field  in  the  vicinity  of  the  transducer. 

Finally,  the  projector  used  to  drive  the  column  into  resonance 
does  not  produce  a  sound  field  identical  to  that  generated  by  a  small 
piezoelectric  transducer  suspended  in  the  water  column.  The  local 
sound  field  near  the  cylinder  is  some  combination  of  radial  modes  which 
matfch  the  cylindrical  boundary  conditions  to  the  overall  plane  wave, 
coupling  the  driving  transducer  to  the  standing  wave  system.  The 
effects  of  this  coupling  factor  are  slightly  frequency-dependent,  and 
are  not  accounted  for  in  this  model. 

Our  measurements  demonstrate  that  the  simple  model  derived, 
despite  these  deficiencies,  produces  consistent  results  to  within  ±0.8 
dB,  but  does  not  explain  differences  from  a  standard  comparison  cali¬ 
bration  of  up  to  2.1  dB.  The  root-mean-square  average  difference  is 
1.5±0.4  dB  for  the  range  of  hydrophone  types  and  frequencies  tested. 
A  method  and  apparatus  for  performing  this  type  of  calibration  experi¬ 
ment,  and  the  results  of  our  experiments,  follow. 


II.  COMPLIANT  TUBE  RECIPROCITY  EXPERIMENTS 


A.  RESONANCE  PARAMETER  MEASUREMENT 

The  formula  for  the  reciprocity  parameter  derived  in  the  preceding 
section  requires  the  accurate  measurement  of  three  resonance  parame¬ 
ters:  the  resonance  frequency,  the  sound  speed  at  that  frequency, 

and  the  resonance  quality  factor:  Each  of  these  measurements  involves 
finding  an  extremum  from  a  set  of  data  samples,  which  will  inevitably 
contain  noise.  The  experimenter  who  attempts  to  resolve  extrema  from  a 
discrete  sample  set  by  comparing  measured  signal  outputs  performs  the 
discrete  equivalent  of  differentiation,  which  enhances  noise  and 
increases  net  error  in  the  result. 

The  solution  to  accurate,  reproducible  resonance  parameter  meas¬ 
urement  which  D.  Conte  and  S.  Garrett  originally  developed  [Ref.  12] 
combines  the  high  data  rate  of  automated  measurements  with  the  power 
of  curve  fitting.  The  intuitive  explanation  of  the  effectiveness  of  this 
method  is  the  human  analogue;  rather  than  take  discrete  measurements, 
experimenters  instead  "swing"  the  system  slowly  through  the 
extremum,  and  let  the  eye  interpolate  the  location  of  that  point.  The 
optimum  automated  processing  analogue  for  this  method  is  curve  fitting 
to  the  minimum  expansion  of  any  function  y(x)  which  will  resolve 
extrema,  the  second  order  equation 


y(x)  =  C0  +  C,x  +  C2x* 


(2.1) 


The  three  coefficients  of  this  curve  can  be  found  for  any  set  of  three 
or  more  measurements;  if  more  than  three  measurements  are  made,  the 
nonlinear  regression  methods  described  by  Bevington  [Ref.  13]  can  be 
used  to  reduce  the  noise  error  in  the  measurement  of  the  extremum,  as 
long  as  the  span  of  the  sample  set  around  the  extremum  is  small  enough 


'  .v  v  .--  - 


that  equation  (2.1)  is  a  valid  approximation.  The  regression  equations 
for  the  three  coefficients  are  algebraic,  and  are  derived  in  Appendix  B 
(equations  (B.12),  (B.13),  and  (B.14)).  The  location  of  the  extremum 
(again  from  Appendix  B)  is  displaced  from  the  coordinate  origin  of  the 
sample  set  by 


—  Xo  —  —  C1/2C2 


(2.2) 


In  pur  application  of  this  method  to  measurement  of  the  resonance  peak 
frequency  ,  the  sample  set  is  composed  of  seven  hydrophone  output 
voltage  amplitudes  measured  at  equal  frequency  intervals;  the  sample 
set  origin  is  shifted  until  the  estimated  peak  frequency  falls  within  one 
interval  of  the  center  frequency  of  the  sample  set.  This  is  the  peak 
frequency  locating  algorithm.  Figure  2.1  is  a  example  of  this  technique; 
the  asterisks  mark  data  points  used  to  establish  the  coefficients  for  the 
curve  and  to  compute  the  location  of  the  maximum,  and  points  marked 
with  circles  are  discarded  from  the  final  fit  because  they  are  more  than 
three  steps  from  the  extremum.  The  sound  speed  is  found  by 
sampling  the  pressure  amplitude  at  regular  intervals  in  depth  from  the 
surface,  and  again  shifting  the  set  used  for  curve  fit  until  the  esti¬ 
mated  maximum  or  minimum  is  within  one  step  of  the  center  of  the  set. 
The  positions  thus  located  occur  at  quarter  wavelength  steps,  and  can 
be  fit  to 


zm  =  m\/4  +  z0 


(2.3) 


to  obtain  the  wavelength;  the  sound  speed  is  then  just  the  product  of 
the  resonance  peak  frequency  and  the  wavelength.  This  is  the  sound 
speed  measurement  algorithm. 

The  curve  fit  in  frequency  also  yields  a  reduced-noise  value  for 
the  peak  amplitude  at  resonance,  which  is  essential  to  accurate 


measurement  of  the  half-power  frequencies  required  for  the  Q  measure¬ 
ment  algorithm  (equation  (1.24)3.  Once  the  half-power  amplitude  is 
known,  the  half-power  frequencies  can  be  found  by  sampling  outward 
from  the  peak  frequency  until  a  pair  of  samples  brackets  this  value  at 
the  minimum  frequency  interval  permitted  by  the  equipment  (in  our 
case,  1  Hz).  The  half  power  frequency  is  then  calculated  by  linear 
interpolation . 

The  curve-fitting  method  can  resolve  extrema  in  noise  to  a  degree 
far  surpassing  the  eye  integration  method  it  was  drawn  from.  It  has 
the  additional  advantages  of  sound  mathematical  basis,  speed,  reprodu¬ 
cibility,  and  independence  from  observer  biases.  It  also  has  the  consid¬ 
erable  advantage  of  being  easier  on  the  observer.  The  method's  main 
drawback  is  its  bias  toward  finding  the  parabolic  curve  in  any  data, 
including  pure  noise;  care  must  be  taken  in  selecting  stepsizes  to 
ensure  that  the  changes  in  the  signal  are  not  so  small  that  they  are 
overwhelmed  by  noise,  and  to  ensure  that  the  sample  set  span  is  not 
larger  than  the  region  about  the  extremum  where  a  second  order 
approximation  is  valid.  The  best  method  for  checking  new  cases  is 
graphical  display  of  the  data  set  and  the  curve  for  the  operator 
approval . 

B.  DESCRIPTION  OF  EQUIPMENT 

1.  Control  System  (Figure  2.3) 

The  experiment  is  controlled  by  a  venerable  Hewlett  Packard 
HP-85  laboratory  computer  with  32K  memory  extension,  IEEE-488  bus 
controller,  and  appropriate  ROM  extensions.  The  program  used  is 
documented  in  Appendix  C.  Commands  generated  by  the  program  over 
the  data  bus  control  the  experiment’s  configuration  changes  during 
measurements  through  the  two  dual  remote  coaxial  switches  (Hewlett 
Packard  HP  59307)  and  a  Hewlett  Packard  HP  3421  Data  Acquisition 
Control  Unit.  The  latter  unit  is  configured  with  two  actuator  relays 
and  eight  measurement  channels  (option  020).  The  two  actuator  relays 
control  the  operation  of  the  reversing  DC  manipulator  drive  motor. 
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Voltmeter 


The  manipulator's  screw  drive,  which  is  described  below,  also 
drives  a  Hewlett  Packard  HED  680  position  encoder  that  increments  a 
Hewlett  Packard  HP  5016  universal  counter  in  totalizing  mode.  The 
changes  in  this  totalized  count  provide  (via  the  data  bus)  the  program 
with  a  measurement  of  actual  manipulator  motion. 

2.  T ransmitter  (Figure  2.4) 

A  Hewlett  Packard  HP  3314  function  generator  provides  a 
stable  bus-controlled  sine  input  signal  at  two  voltage/current  levels, 
and  a  synchronization  signal  at  the  selected  frequency  and  constant 
volt'age.  The  high  voltage  output  is  amplified  and  used  to  drive  the 
piezoceramic  transducer  at  rms  voltages  of  10  to  20  volts  AC.  The 
low-voltage,  high-current  demands  of  the  USRD  Jll  projector  [Ref.  14] 
are  met  by  low-gain  power  amplification  of  the  function  generator's  low 
voltage  output.  The  transmit  remote  coaxial  switch  (714) 2  selects  the 
signal  source  from  these  two  options  and  feeds  the  signal  to  the 
reversible  transducer  switch  (706A),  to  the  Jll  projector,  or  through  a 
test  attenuator  which  returns  a  known  signal  to  the  receiver  input 
switch  (706B). 

A  Hewlett  Packard  HP3478  digital  multimeter,  configured  as  a 
remotely  controlled  5  1/2  digit  AG  voltmeter,  is  connected  to  the 
transmit  signal  line  to  measure  the  transmitting  voltage. 

3.  Receiver  (Figure  2.4) 

The  test  transducers  used  in  these  experiments  are:  two 
Celesco  LC-10  piezoceramic  hydrophones  (Figure  2.5),  selected  for  their 
small  size  and  good  transmitting  response;  and  a  Celesco  LC-32  piezo¬ 
ceramic  hydrophone  (Figure  2.6)  of  larger  size.  The  voltage  sensitivi¬ 
ties  of  all  three  units  were  measured  by  comparison  with  a  reference 
standard  USRD  H56  hydrophone  [Ref.  14].  The  results  are  listed  in 
Table  II.  The  LC-10's  are  mounted  on  0.375  inch  thinwall  copper 
tubing;  the  LC-32  is  mounted  on  0.4375  inch  tubing.  Signal  cables  are 


2The  remote  switches  are  identified  here  by  their  IEEE-488  bus 
address  codes. 


TABLE  II 

Comparison  Calibration  Results 
(All  values  £  0.8  decibels) 


Frequency 

LC-10#2166 
&  2m  cable 

LC-10#2338 
&  9m  cable 

LC-32#85 
&  11m  cable 

740 

-206.5 

-206.8 

-212.3 

860 

-206.8 

-207.3 

-213.4 

1000 

-206.9 

-208.8 

-214.3 

routed  up  the  interior  of  the  tubing  to  minimize  effects  on  the  sound 
field.  The  two  test  transducers  are  suspended  from  the  manipulator 
platform.  Their  output  is  returned  from  the  top  of  the  manipulator 
assembly  by  two  of  three  coaxial  cables  of  the  same  length  and  meas¬ 
ured  capacitance.  The  third  of  these  cables  connects  the  test  attenuator 
into  the  receiver  input  switch  (706B),  which  makes  all  receive  signal 
path  wiring  nearly  equivalent.  The  reversible  transducer  switch  (706A) 
selects  either  its  receive  or  transmit  lines;  the  receiver  input  switch 
(706B)  selects  the  lockin  analyzer  input  from  the  reversible  transducer 
receive  line,  the  test  hydrophone,  or  the  test  attenuator. 

The  EG&G  PAR  model  5204  lockin  analyzer  filters  out  the 
components  of  the  received  signal  which  are  in-phase  and  in-quadrature 
with  the  reference  (synchronization)  signal.'  Their  Pythagorean  sum  is 
used  to  obtain  the  magnitude.  The  analyzer  is  typically  operated  with 
an  integration  time  of  300  milliseconds,  which  corresponds  to  an  equiva¬ 
lent  noise  bandwidth  of  0.8  Hz;  this  causes  the  lockin  analyzer  to  act 
as  a  very  sharp  filter  synchronized  to  the  transmitter  signal.  The 
analog  output  of  the  model  5204  is  passed  to  a  measurement  channel  of 
the  HP  3421  for  digitization. 

4.  Mechanical  System 

The  tests  are  conducted  in  a  6-inch  (154  mm)  inside  diameter 
Schedule  40  polyvinyl  chloride  tube  mated  to  the  active  surface  of  a 
USRD  J11  electrodynamic  underwater  projector  [Ref.  14].  The  tube 
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DIMENSIONS  IN  INCHES 


Figure  2.5  Celesco  LC-10  Hydrophone 
(Reprintecf  with  permission  of  Celesco  Industries, 


Inc. ) 


stands  1180  mm  above  the  ill’s  active  surface,  and  is  filled  to  a  depth 
of  about  1  meter;  its  wall,  which  is  smooth  on  the  interior  except  for 
small  regular  extrusion  marks,  is  8  mm  thick.  The  configuration 
resembles  the  USRD  G40  calibrator  used  for  comparison  calibrations 
[Ref.  14].  (This  system  was  originally  constructed  by  Mills  and 
Garrett  for  their  work  [Ref.  15]  with  fiber  optic  gradient  hydro¬ 
phones.)  Equalization  for  the  J11  driver  is  provided  by  air  pressuriza¬ 
tion  of  the  installed  air  bag,  monitored  with  a  water-filled  U-tube 
manometer. 

The  test  transducers  are  suspended  on  thinwall  tubing 
supported  on  a  rigid  table  which  is  raised  and  lowered  on  two  one-inch 
rod  tracks  by  a  motor-driven  screw  drive;  maximum  range  of  motion  is 
800  mm.  The  drive  controls  were  described  in  the  paragraphs  on  the 
control  system  above. 
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Figure  2.6  Celesco  LC-32  Hydrophone 
(Reprinted  with  permission  of  Celesco  Industries,  Inc.) 


C.  PROCEDURES  FOR  EXPERIMENT 


1.  Preparation 


Before  the  test  transducers  are.  suspended  in  the  water 
column,  the  experimenter  must  measure  the  reversible  transducer's 
input  capacitance,  so  that  the  transducer  current  can  be  computed  from 
the  transmitter  input  voltage  and  drive  frequency.  The  active  surface 
lengths  of  both  hydrophones  should  also  be  measured  for  use  in 
computing  the  transducer  amplitude  averaging  factor  (equation  (1.33)). 
The  test  transducers  are  then  suspended  from  the  manipulator,  wetted 
with  a  detergent  solution,  immersed  in  the  water  column,  and  brushed 
to  remove  bubbles  from  the  active  surface.  If  necessary,  the  tube's 
inner  wall  should  also  be  brushed  free  of  bubbles.  The  water  column 
depth  is  measured  with  a  meter  stick. 


2.  Mode  Selection 


Two  parameters  are  available  to  vary  the  resonant  frequency 
of  the  water  column  at  which  the  calibration  is  done:  the  water  depth 
L  and  the  mode  number  n  (equation  (2.4)).  The  procedure  requires 
wavelengths  much  larger  than  the  coaxial  lengths  of  the  transducers. 
In  particular,  the  wavelength  should  be  much  longer  than  the  coaxial 
length  of  the  reversible  transducer  T  to  ensure  that  T's  surface  is 
entirely  within  the  velocity  node  and  feels  little  or  no  impressed  pres¬ 
sure  when  T  is  transmitting.  The  wavelength  should  also  be  at  least 
four  times  the  length  of  the  test  hydrophone  H.  Practical  experience 
has  shown  that  calibration  results  are  most  accurate  when  T  and  H  are 
positioned  at  two  different  pressure  maxima  of  the  n=3  or  higher  order 
resonance  modes. 


L  =  |(n  +  1/2)  =jf  (n  +  1/2) 


(2.4) 


3.  Resonance  Parameter  Measurements 

After  positioning  the  test  transducers  at  the  approximate  loca¬ 
tions  of  the  deepest  and  next  deepest  accessible  spatial  pressure  maxima 
with  similar  pressure  amplitudes,  the  operator  uses  the  peak  frequency 
measurement  algorithm  to  find  the  resonance  frequency  f.  The  trans¬ 
ducers  are  then  withdrawn  above  the  first  pressure  minimum. 

One  of  the  transducers  is  now  used  to  measure  the  locations 
of  maxima  and  minima  in  the  standing  wave  pattern.  Linear  regression 
of  this  data  gives  the  quarter  wavelength,  which  in  turn  gives  the 
sound  speed. 

Finally,  the  transducers  are  precisely  repositioned  at  the 
maxima  where  the  resonance  frequency  measurement  was  made.  A 
precise  measurement  of  the  peak  frequency  and  of  the  Q  factor  are 
made,  completing  the  set  of  measurements  needed  to  compute  the  reci¬ 
procity  parameter  from  equation  (1.32). 


4.  Calibration  Measurements 

Before  making  the  final  measurements,  the  experimenter 
switches  the  attenuated  test  signal  into  the  receiver  output  to  measure 
the  net  effects  of  cable  losses,  preamplifier  gain  drift,  and  deviations 
between  the  lockin  analyzer's  and  the  digital  voltmeter's  measurements. 
This  correction  is  applied  along  with  the  amplitude  averaging  correction 
factors  from  equation  (1.33)  to  the  calibration  voltage  measurements 
that  follow. 

The  transmitter  is  switched  to  the  reversible  transducer  T, 
and  the  receiver  to  the  test  hydrophone  H.  The  transducer-to- 
hydrophone  output  voltage  E(TH)  is  measured,  corrected,  and 
recorded,  as  is  the  transmit  voltage  V(T).  The  latter  is  used  with  the 
capacitive  admittance  of  T  to  compute  the  transmitter  current  l(T). 

The  transmitter  is  then  switched  to  the  projector,  and  the 
drive  voltage  adjusted  to  give  a  sound  field  of  the  same  amplitude  at  H 
as  was  used  in  the  previous  step.  The  projector-to-hydrophone  output 
voltage  is  recorded  (it  is  equal  to  E(TH),  within  the  limits  of  the 
projector  drive  adjustment  error).  The  final  calibration  measurement  is 
made  with  the  receiver  switched  to  T  with  the  projector  drive 
unchanged;  the  result  is  the  projector-to-transducer  output  voltage 
E(PT) . 

These  measurements  complete  the  requirements  to  compute  the 
test  hydrophone  voltage  sensitivity  from  equation  (1.5).  The  same  set 
of  measurements  can  be  used  to  compute  the  reversible  transducer's 
sensitivity  using 


Mt  = 

Eph’t 


(2.5) 


The  results  for  M  are  in  volts  per  pascal;  the  customary  units  quoted 
are  decibels  re  1V/yPa,  which  is  found  by  calculating 


Jt  =  20  log  M  -  120 


(2.6) 
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5.  Self-consistency  Check 


When  both  transducer  can  be  used  as  reversible  transducers, 
the  results  of  the  reciprocity  calibration  can  be  checked  promptly  by 
exchanging  the  transducers'  roles  in  the  calibration  measurements. 
Repeat  measurement  of  the  resonance  parameters  and  transducer  reposi¬ 
tioning  is  not  required. 

Unless  the  two  transducers  have  identical  nonlinear  responses, 
this  second  set  of  readings  will  confirm  the  accuracy  of  the  assumption 
of  reciprocity.  It  will  not,  however,  measure  the  error  in  computing 
the.  reciprocity  parameter  J,  which  must  be  estimated  independently. 

D.  RESULTS  OF  EXPERIMENTS 

1 .  Reciprocity  Tests 

Tests  were  run  on  the  three  transducers  in  all  pair  combina¬ 
tions  at  the  three  accessible  modes:  n=  3,  4,  and  5.  Modes  below  n=3 
were  inaccessible  because  the  LC-10  transducers'  transmitting  response 
was  too  small  for  reliable  measurements;  modes  above  n=5  were  attenu¬ 
ated  into  the  ’tails’  of  the  lower  resonances,  which  obscured  the  half¬ 
power  points  needed  to  measure  Q.  '  Comparison  of  a  number  of 
reciprocity  calibration  test  results  showed  reproducible  results  to  within 
0.8  decibels.  However,  the  calibration  results  in  Table  IV  disagree 
systematically  with  the  values  obtained  for  the  same  hydrophones  by 
comparison  with  the  USRD  standard  hydrophone  (Table  II)  by  1.5  deci¬ 
bels  (Table  V) . 

2.  Error  Analysis 

The  self-consistency  error  quoted  in  Table  IV  was  determined 
by  reversing  the  electrical  roles  of  the  two  transducers  without 
changing  their  positions;  the  difference  in  computed  results  for  the  two 
measurements  is  the  degree  to  which  the  system  is  not  reciprocal.  The 
estimated  error  is  the  combination  of  the  self-consistency  error  with  the 
error  in  measuring  the  reciprocity  parameter,  listed  in  Table  III. 


The  primary  source  of  self-consistency  error  is  the  difference 
in  pressure  amplitudes  at  the  two  maximum  locations  used,  which  we 
believe  is  caused  by  attenuation.  By  careful  selection  of  modes  and 
locations,  we  were  able  to  reduce  this  error  to  0.2  dB.  Another  lesser 
source  of  self-consistency  error  is  the  changing  effect  of  the  support 
tubing  used  for  the  transducers  on  the  sound  field  as  they  are  lowered 
into  the  water  column.  Insertion  of  the  transducers  and  their  supports 
into  the  water  column  displaces  about  0.4  percent  of  the  water  column 
in  volume,  causing  a  similar  change  in  the  column  depth  and  in  the 
positions  of  maxima  and  minima.  The  change  in  the  minima  is  the  most 
significant  error  encountered  in  measuring  the  wavelength  of  the 
standing  wave;  because  the  cosine  function  is  flat  to  first  order, 
however,  the  effect  on  transducer  positioning  at  maxima  is  negligibly 
small. 

Compliant  materials  associated  with  the  signal  cables  and  the 
transducers  appear  to  degrade  the  reciprocal  nature  of  the  geometry. 
This  error  is  reduced  to  very  low  values  by  selection  of  transducer 
location  and  by  eliminating  compliant  materials  (signal  cables,  for 
example)  in  contact  with  the  water. 

Parameter  measurement  error  is  dominated  by  the  apparent 
variation  in  the  system  quality  factor  over  the  length  of  the  tube.  An 
increase  of  6-10  percent  in  the  value  of  Q  determined  by  the  half-power 
frequency  method  was  observed  as  the  depth  of  the  measurement  was^a., 
increased  by  one  wavelength;  our  model  for  the  compliant  tube  system 
does  not  account  for  such  effects,  and  is  therefore  limited  in  the  accu¬ 
racy  with  which  it  predicts  the  reciprocity  parameter  to  about  0.8  dB. 
This  is  the  largest  recognized  source  of  measurement  error  in  the 
compliant  tube  reciprocity  method. 

The  systematically  lower  sensitivities  measured  by  the  reci¬ 
procity  method  relative  to  those  measured  by  comparison  with  the  H56 
standard  are  not  yet  understood.  This  discrepancy  may  rise  from  the 
inadequacies  of  our  model  for  the  system  energy  balance,  from  the 
method  used  to  measure  Q,  or  even  from  technical  difficulties  with  the 
comparison  calibration. 


Measurement  of  Wall  Motion 


Since  sound  speed  can  be  measured  directly,  it  is  not  neces¬ 
sary  to  attempt  the  difficult  measurement  of  the  actual  elastic  strain  in 
the  tube  wail  as  part  of  the  calibration  process.  However,  the  degree 
of  correspondence  of  the  actual  strain  behavior  to  the  model  derived  in 
the  first  chapter  is  an  important  indication  of  the  validity  of  our 
method . 

The  circumferential  wall  strain  was  measured  by  attaching 
constantan  .250  inch  strain  gages  [Ref.  16]  to  the  exterior  surface  of 
the,  tube  at  the  depths  of  some  of  the  pressure  maxima  located  by  the 
test  transducer.  A  second  lockin  analyzer  used  as  the  sensing  element 
of  a  resistive  DC  bridge  measured  the  magnitude  and  phase  of  the 
strain  oscillations  sensed  by  the  strain  gage;  by  comparing  the  output 
of  the  strain  gage  with  the  magnitude  and  phase  of  the  pressure  meas¬ 
ured  by  the  transducer  at  that  depth,  we  computed  the  complex  impe¬ 
dance  (W,9)  of  the  wall  ring  (equation  (1 . 21 ))  .  The  results  in  Table  VI 
agree  in  general  with  the  measured  sound  speed  in  the  tube.  Also,  the 
indicated  damping  phase  effects  are  reasonably  small,  justifying  our 
approximation  for  the  wall  potential  energy  computation  of  equation 
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Variations  in  Calibration  Results 
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The  upper  frequency  limit  depends  both  on  the  onset  of  azimuthal 
modes,  and  on  the  increasing  attenuation  of  the  sound  wave  by  the  PVC 
wall.  The  former  constraint  breaks  the  plane  wave  symmetry  when  the 
frequency  is  greater  than  (0.586*L/R)  times  the  fundamental  longitu¬ 
dinal  mode  [Ref.  17].  The  attenuation  problem  causes  excessive  varia¬ 
tion  of  the  measureable  system  quality  factor  with  depth  and  reduces  Q 
to  the  point  where  half  power  frequencies  are  obscured  by  adjacent 
resonances . 

Test  transducer  geometry  is  restricted  to  shapes  that  are  small 
enough  in  cross  section  that  reflections  from  their  lower  surface  do  not 
significantly  affect  the  standing  wave.  The  length  of  the  reversible 
transducer,  as  discussed  in  the  procedure  description,  should  be  as 
short  as  possible  consistent  with  adequate  transmitting  response.  (A 
flat  disk  transducer  might  be  more  effective  than  the  cylindrical  LC-10 
we  used.)  The  test  hydrophone  can  be  much  larger,  if  the  transducer 
amplitude  averaging  correction  is  properly  computed. 

Accuracy  of  the  method  is  primarily  limited  by  the  effects  of  the 
wall  attenuation.  A  derivation  which  accounted  for  these  effects,  and 


for  the  lesser  effects  of  disturbances  to  the  standing  wave  caused  by 
the  test  transducers  and  their  support  tubing,  could  conceivably 
achieve  much  lower  calibration  uncertainities . 

A  significant  advantage  of  the  compliant  tube  method  is  the  use  of 
pressure  antinodes  as  virtual  rigid  boundaries  at  which  to  locate  the 
test  transducers.  Many  of  the  difficulties  associated  with  mounting  a 
transducer  in  an  real  boundary  which  must  have  some  compliance  are 
thus  avoided. 

B.  POTENTIAL  APPLICATIONS 

Two  applications  for  the  compliant  tube  system  on  a  practical  scale 
have  been  suggested  to  the  author:  a  field/small  lab  calibrator,  and  a 
large,  very  low  frequency  calibrator. 

The  portability  and  flexibility  of  the  test  system  used  for 
compliant  reciprocity  recommends  the  method  for  general  calibration  use 
in  small  laboratories  and  in  the  field.  The  1 -meter  PVC  tube  provides 
calibration  capability  in  the  1  kHz  range;  a  2-meter,  10-inch  diameter 
tube  would  be  more  versatile  and  could  reach  frequencies  of  500  Hz. 
The  manipulator  assembly  need  not  be  as  complicated  or  massive  as  the 
one  used  in  our  experiment;  any  system  which  can  position  the  trans¬ 
ducers  within  a  few  millimeters  is  adequate.  Because  the  method 
involves  resonance  measurement  methods  as'  well  as  the  techniques  of 
reciprocity  calibration,  it  may  also  be  of  interest  for  advanced  instruc¬ 
tional  laboratories. 

On  a  larger  scale,  a  10-meter,  30-inch  diameter  tube,  which  would 
fit  upright  in  a  two-story  open  bay,  could  reach  frequencies  below  100 
Hz  and  calibrate  units  similar  in  size  to  the  Navy's  TR-155  transducer. 
If  current  interest  in  very  low  frequency  underwater  acoustics  indicates 
a  need  for  an  alternative  to  the  traveling-wave  tube  system  described 
in  [Ref.  8],  compliant  tube  reciprocity  may  provide  the  answer  at  a 
very  low  cost. 

There  is  also  no  reason  we  are  aware  of  that  the  benefits  of 
compliant  waveguide  usage  should  be  limited  to  resonant  calibration. 


Construction  of  a  compliant  traveling  wave  tube  calibrator  with  a  non¬ 
reflecting  termination  or  of  a  Helmholtz  resonant  calibrator  using  two 
compliant  tubes  should  also  be  practical.  The  Helmholtz  calibrator 
concept  in  particular  has  the  potential  of  extending  the  range  of  accu¬ 
rate  absolute  calibration  for  the  underwater  acoustician  to  extremely  low 
frequencies. 


APPENDIX  A 
GLOSSARY 


Subscripted  variables  in  the  text  are  represented  using  paren 
theses,  due  to  word  processor  limitations. 

A  Cross  sectional  area  of  the  water  column  at  rest. 

A(T)  Active  surface  area  of  the  transducer, 

a  Instantaneous  cross  sectional  area  of  the  water  column. 

C(0  Parabolic  curve  fit  coefficients, 

c  Sound  speed  in  a  compliant  tube. 

D  Distensibility  of  the  boundary  wall. 

£(TH)  Transducer-driven,  test  hydrophone  output  voltage. 

E(PH)  Projector-driven,  test  hydrophone  output  voltage. 

E(PT)  Projector-driven,  reversible  transducer  output  voltage, 

f  Resonance  frequency. 

H  The  test  hydrophone  used  only  to  receive  during  the  reci 

procity  calibration. 

h  Thickness  of  the  tube  wall. 

I  (T)  Transmitter  current  in  reversible  transducer  T. 

J  Reciprocity  parameter,  or  acoustic  transfer  admittance 

between  two  points. 

K  Adiabatic  compressibility  of  water, 

k  Wave  number  at  resonance  =  2it/X. 

L  Depth  of  the  water  column. 

L(T)  Transducer's  active  length  normal  to  the  plane  wave. 

M  Hydrophone  voltage  sensitivity. 

m  Number  of  quarter  wavelengths  below  the  water  surface, 

n  Resonance  mode  index. 

n  Unit  vector  normal  to  surface  element  ds. 

p  Acoustic  wave  excess  pressure. 

P  Pressure  amplitude  of  the  acoustic  wave. 

Q  Resonance  quality  factor. 

R  Inner  radius  of  the  tube. 


i 


3 


R(m) 


x(0) 
x  (m) 
Y 

y(i) 

Z(E) 

Z(m) 


Damping  coefficient  for  tube  wall  vibrations. 

Transmitting  current  response. 

Surface  area  element. 

The  test  transducer  used  both  as  a  projector  and  as  a 
hydrophone  during  the  reciprocity  calibration. 

Volume  velocity  of  the  acoustic  wave. 

Particle  velocity  of  the  acoustic  wave. 

Volume  element. 

Complex  mechanical  impedance  magnitude  of  the  tube  wall. 
Center  coordinate  of  the  data  set  used  for  curve  fitting. 
Coordinate  of  extremum  of  the  curve. 

Elastic  coefficient  for  tube  wall  expansion. 

Sample  set  to  be  curve  fit. 

Transducer  blocked  electrical  impedance. 

Transducer  open-circuit  mechanical  impedance. 

Depth  from  surface  of  water  column. 

Depth  measurement  offset. 

Depth  of  the  mth  pressure  maximum/minimum. 

Sampling  interval  for  curve  fit  sample  set  y(i). 

Areal  strain  in  the  wall.  See  equation  (1.19). 

Length  averaging  correction  factor. 

Complex  mechanical  impedance  phase  angle  of  the  tube  wall. 
Wavelength  of  the  standing  wave. 

Instantaneous  local  water  density. 

Ambient  water  density. 

Transformation  factor. 

Sum  of  squared  errors  in  a  curve  fit  to  data. 

Angular  frequency  at  resonance. 

Areal  expansion  coordinate. 


GSWP 


APPENDIX  B 

PARABOLIC  SMOOTHING  ALGORITHM 


The  algebra  of  deriving  formulas  for  the  best  choice  of  coefficients  of  a 
nonlinear  equation,  such  that  the  sum  of  squared  errors  described  by 

xl  =  Exf=E 

is  minimized  ,  is  described  by  Bevington  in  [Ref.  13].  The  minimum 
for  the  function  defined  by  equation  (B.1)  is  defined  by  the  selection 
of  coefficients  which  solve 


i 

y(xt)  -  y« 


(B.1) 


ax! 

ac 


=  o 


( B .  2) 


where  C  is  any  of  the  coefficients  of  the  curve  to  be  fit  to  the  data. 

The  simplest  function  that  will  predict  maxima  and  minima  is  the 
quadratic  equation 


y(x)  =  Co  +  CiX  +  C2X2 


(B  .3) 


The  single  extremum  of  this  curve  occurs  at 


xm  x  0  —  C  i  /2  C  2 


( B .  4) 


with  respect  to  the  coordinate  origin  located  at  x.  Let  the  sample  set 
be  seven  measurements  of  the  function  y  sampled  at  regular  intervals  5 


yi  *  y(x0  +  i<5)  ;  i=  -3,  -2.  -1,  0,  +1,  +2,  +3. 


(B.5) 


Evaluating  the  best  fit  function  s  partial  derivatives  with  respect  to  the 
three  parameters  for  such  a  sampled  data  set  gives  equations  (B.6) 
thru  (B.11);  summing  and  cancelling  yield  equations  (B.12),  (B.13), 
and  (B.14).  These  analytic  results  can  then  be  used  with  equations 
(B.3)  and  (B.4)  to  obtain  the  location  of  the  maximum  and  its  ampli¬ 
tude,  or  the  smoothed  value  for  the  amplitude  at  any  location. 


X-3  —  95rC2—  35C,+  Co-  y-3 

(B.6) 

X-2  ~  45  C 2  —  25C,+  Co-  y -2 

(B  .7) 

X-i  —  (52C2-  5Ci+  Co-  y -1 

(B  .8) 

X  0  =  Co~  y  0 

(B.9) 

X+\=  <5zC2+  Co- y+( 

X+2  =  45  C2+  25Ct+  Cq—  y+z 

(B. 10) 

X+i  =  95  C2  +  35C,+  Cq—  y+3 

(B.ll) 

(B.12) 

(B.13) 


84<52C2  =  5(y*j+  y_3)  -  3(y<.,+  y_,-  4y0  ) 

28<5C,  =  3(y+i-  y.})  +  2(y+2-  y-2)  +  (y+,-  y_,) 

21C0  =  7y0  +  6(y*i  +  y_,)  +  3(y+2+y_2)  -  2(y+J+  y_3) 


(B.14) 


APPENDIX  C 

ACTRE:  HP-85  BASIC  PROGRAM  FOR  RECIPROCITY  EXPERIMENTS 


A.  WHAT  THE  PROGRAM  DOES 

The  Automated,  Compliant  Tube  Reciprocity  Experiment  (ACTRE) 
program  provides  the  operator  with  a  menu  of  soft  key  controlled 
"tools"  for  resonance  parameter  and  calibration  measurements.  Its 
des.ign  was  based  on  the  need  for  maximum  flexibility,  rather  than 
maximum  unattended  operation. 

All  of  the  measurement  "tools"  are  subroutines  called  from  the 
HP-85  function  keys  by  the  operator.  The  major  procedures  are:  the 
peak  frequency  location  algorithm  (KEY  1);  the  spatial  optimum  location 
algorithm,  which  also  computes  wavelength  and  sound  speed  (KEY  2); 
the  Q  measurement  algorithm  (KEY  5);  and  the  reciprocity  calibration 
measurement  sequence  (KEY  7).  Supporting  procedures  provide  the 
actual  control  statements  for  remote  switching,  function  generator 
setup,  manipulator  positioning  of  the  hydrophones,  receiver  output 
measurement,  and  graphics  output.  Any  procedure  may  be  terminated 
by  calling  the  process  clear  sequence  (KEY  8).  Comments  on  the 
detailed  structure  of  the  program  are  included  with  the  source  code. 

The  program  was  developed  to  fit  the  demands  of  the  experimental 
apparatus, the  experimenter,  and  the  research  goals  of  this  work  as  well' 
as  the  limits  of  the  host  computer  permitted.  Time  and  core  memory 
constraints  severely  limited  efforts  to  make  the  program  friendly  or 
accessible  to  the  uninformed. 

B.  TEST  EQUIPMENT  BUS  CONNECTIONS 

The  code  is  designed  to  use  the  IEEE-488  bus-controlled  test 
equipment  described  below.  All  of  the  primary  measurements  are 
made  with  bus-controlled  equipment,  except  for  the  hydrophone  output 
voltage  measurement  by  the  EG&G  5204  lockin  analyzer,  whose  output 
signal  must  be  converted  using  the  HP  3421  D/A  unit. 


Bus  assignment  numbers  used  in  the  program: 


705  HP7470  hardcopy  plotter. 

706  HP59307A  remote  controlled  switch: 

A  switch:  reversible  transducer; 

B  switch:  receiver  input. 

707  HP  3314  function  generator. 

709  HP  3421  Data  Acquisition/Control  Unit,  with  Opt  020 
board  in  slot  0.  Channel  assignments: 

0:  Manipulator  motor  drive  up  actuator 
1 :  Manipulator  motor  drive  down  actuator 
5:  5204  lockin  analyzer  magnitude  signal. 

714  HP  59307A  remote  controlled  switch  (transmit  side). 

720  HP  5316  Frequency  counter,  in  count  totalize  mode. 

723  HP  3478  digital  multimeter,  configured  as  an  AC  volt¬ 

meter. 


C.  VARIABLE  ASSIGNMENTS 
1 .  Defined  Functions 

FNH(XI)  Length  averaging  correction  function  (equation 

(1.33)). 

FNK(KO)  Cycles  K0  through  the  indices  1,2,3.  Used  in  plot 
scaling  algorithm. 

FNL(XO)  Returns  decibel  value  of  X0  =  20  log(XO). 

FNM(XO)  Returns  voltage  sensitivity  in  dB  re  1  V/uPa  (equa¬ 
tion  (2. 6)  ) . 

FNR(XO)  Returns  value  of  X0  rounded  to  three  significant 
places;  required  for  compatibility  with  HP3314  function 
generator  amplitude  command. 


2.  Variables 


A(7) 

A2 


A4 

A5 


A7 

A8 


Amplitude  sampling  array,  used  for  curve  fitting  algo¬ 
rithms  . 

Amplitude  ratio  of  hydrophone  output  voltage  to 
transmitter  input  voltage.  Measured  by  the  Amplitude 
Read  procedure. 

Half  power  amplitude  in  Q  measurement  procedure. 
Lockin  analyzer  output  signal  to  D/A  unit  (709) 
channel  5;  ranges  from  0-1  V  to  indicate  0-full  scale 
magnitude  reading. 

Resonance  peak  amplitude,  computed  from  curve  fit. 
Preceding  measurement  storage--in  Q  measurement 
procedure. 


B(3)  Plot  scaling  parameter  array 

B1-B9  Plot  parameters 


C1-C5  Sums  used  in  wavelength  measurement  least-squares 
procedure  [Ref.  13]. 

Cl :  Sum  of  N9. 

C2:  Sum  of  N9*N9. 

C3:  Sum  of  L9. 

C4:  Sum  of  L9*N9. 

C5:  Sum  of  L9*L9. 

C6  The  quantity  A  in  the  least-squares  fit. 

C7  Resonance  sound  speed  (m/s). 

C3  Error  estimate  for  C7  (m/s). 


D1 

Water  column 

depth(mm) 

D2 

Water  column 

cross  section  (mm2) 

D3 

Water  column 

volume(m3 ) 

D5 

Length  of  T 

(mm) . 

D6 

Length  of  H 

(mm) . 
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El  E(PT):  Output  voltage  of  T  ,  P  driven  (V). 

E2  E ( PH ) :  Output  voltage  of  H  ,  P  driven  (V). 

E3  £(TH):  Output  voltage  of  H  ,  T  driven  (V). 

All  frequencies  are  in  hertz(Hz). 

FO  Center  frequency  of  sample  set. 

FI  First  peak  detection  /  Lower  1/2  power  point. 

F2  Second  peak  detection  /  Upper  1/2  power  point. 

F3  Peak  frequency  estimate. 

F4  Lower  search  boundary. 

F5  Upper  search  boundary. 

F6  Drive  frequency;  argument  passed  to  Amplitude  Read 

Setup  procedure. 

F7  Resonance  peak  frequency. 

F8  Error  estimate  for  F7. 

F9  Stepsize. 

G2  Receiver  preamplifier  gain. 

G5  Lockin  analyzer  output  expand  setting. 

G7  Calibration  attenuator  value.' 

G8  Drive  boost  factor  for  preamp  calibration  test  in 

Calibration  routine. 


HI 

H2 


Length  averaging  correction  factor  for  T. 
Length  averaging  correction  factor  for  H. 


General  use 

Size  of  plot  files  X(),  Y(). 

Selection  index  for  Operator  Input  procedure. 

Number  of  time  constants  T5  to  delay  before  measure¬ 
ments  . 

Select  index  for  Switching  procedure. 

Wavelength  measurements  count. 

Reversible  transducer(T)  transmit  current  (A). 
Reciprocity  factor  (V*A/Pa2). 

Water  density  (kg/m5). 

Plot  scaling  array  index. 

Pressure  maximum/minimum  position  (m) . 

Hydrophone  voltage  sensitivity  of  T  (V/Pa). 
Hydrophone  voltage  sensitivity  of  H  (V/Pa). 

Position  counts;  used  to  corrtpute  manipulator  motion. 
Number  of  1/4  wavelengths  from  surface  corresponding 
to  L9. 

Measured  quality  factor  of  resonance. 

Position  encoder  count  ratio  for  down  motion 
(counts/mm) . 

Manipulator  drive  rate  (counts/sec),  down  motion. 
Manipulator  drive  rate  (mm/sec),  down  motion. 

Allowed  projector  drive  adjust  error. 


R4  Allowed  peak  frequency  detection  range. 

R5  Allowed  min/max  detection  range. 

R7  Q  measurement  coarse  stepsize  factor. 

R8  Manipulator  up/down  motion  rate  ratio. 

The  S  variables  are  all  control  flags/indices. 

50  General  use. 

51  Read  enable  flag. 

52  Manipulator  motion  stop  flag. 

53  Plot  option  enable  flag. 

54  Fast  display  refresh  flag. 

55  Drive  voltage  select  index: 

0:  use  J11  projector  drive  voltage  V4; 

1:  use  piezoceramic  drive  voltage  V3. 

56  Automatic  drive  voltage  control  enable  flag. 

57  Resonance  peak  detection  index. 

S9  Plot  execute  flag  and  abscissa  select  index. 

All  times  are  in  milliseconds. 

T2  Manipulator  motion  time  adjust. 

T5  Lockin  analyzer  time  constant. 

T7  Additional  delay  time  before  measurement. 

T8  Size  of  additional  delay  T7  to  be  used  when  required 


U0  Curve  fit  amplitude  at  center  of  sample  set. 

U1  Curve  fit  slope*stepsize. 

U2  Curve  fit  curvature  *(stepsize) 2 . 

U9  Stepsize  used  in  last  curve  fit  (  for  plot  procedure). 


All  voltages  are  in  volts(V)rms,  except  where  indicated. 

Transmitter  input  voltage. 

Hydrophone  output  voltage. 

Piezoceramic  drive  voltage(Vp-p) . 

J11  projector  drive  voltage(Vp-p) . 

Lockin  analyzer  fullscale  setting;  input  as  milli¬ 
volts,  and  converted  to  volts. 

Function  generator  drive  voltage  (Vp-p). 

Transmitter  input  voltage  at  resonance. 

Curve  fit  amplitude  at  F7,  returned  by  5-point 
Voltage  Measurement  procedure. 

W0-W3  General  use  working  variables. 

X(50)  Abscissa  plot  storage  array  (Hz  or  mi*). 

Y(50)  Ordinate  plot  storage  array  (mV  rms). 

XO  Reference  x-coordinate  for  curve  fit  plot  (Hz  or  mm). 

All  variables  beginning  with  "Z"  are  in  millimeters(mm) . 

Z(7)  Sampling  positions  for  X  Measurement  procedure. 

ZO  Manipulator  table  position. 

Z2  Estimated  position  of  pressure  minimum/maximum. 

Z4  Position  of  center  of  sample  data  set/  operator-ordered 

position  to  move  to. 

Z8  Measured  stepsize  in  X  Measurement  procedure. 

Z9  Ordered  stepsize. 

Character  String  Variables . 

A$(128)  Title;  used  to  identify  printout  and  plots 

M$(20)  Code  identifier. 

SS(5)  Remote  switch  position  identifier. 

W$(32)  General  character  input. 


VI 

V2 

V3 

V4 

V5 

V6 

V7 

V8 


D.  SOURCE  CODE 


1000  !  program  actre-7 

1010  !  Automated,  Compliant  Tube  Reciprocity  Experiment 
1020  !  M.8.  Johnson,  Dec  34 
1030  ! 

1040  PRINTER  IS  2  9  OPTION  BASE  19  DES 

1050  DIM  AtC128],H*C5],W*C323  ,SU53 

1060  A$="ACTRE-7  mod  3  "  9  CLEAR  9  DISP  A$ 

1070  i  Functions 

1080  DEF  FNH(Xl)  =  X 1  /S  IN ( RTD  ( X 1 ) ) 

1090  DEF  FNK(Kl)  =  K1  MOD  3+1 
1100  DEF  FNL(Xl)  *  20*LGT(X1) 

1110  DEF  FNM(Xl)  =  FNL ( X 1 ) - 1 20 

1120  DEF  FNRtXl!  1  rounds  to  3  significant  figures  for  HPT 

1130  W l *.01 

1140  !  begin  1 oop 

1150  IF  I P ( X 1 >  <  1  THEN  X 1  =  X  1*1000  0  W1 =W1 / 1000  9  GOTO  1140 
1160  IF  I P ( X 1 ) >9  THEN  X  1  =  X 1  /  1 0  9  W1=W1*10  8  GOTO  1140 
1170  !  end  loop 
1180  FNR*N 1  +  IP ( X l* 100+.  5) 

1190  FN  END 

1200  !  end:  functions 

1210!  Variable  initialization 

1220  REAL  A(7) ,A2,A4,A5,A7,A8,C1,C2,C3,C4,C5,C6,C7,C3ID1 ,D 
1230  REAL  D3,El,E2,E3,F0,Fl,F2,F3,F7,S2,J3,J4,J7,K0,i.9,Ml,l 
1240  REAL  NO , N 1 , N2 , UO , U1 , U2 , V 1 , V2, VB , Y3 , Z ! 7! , 20 , Z2 , 23 , 24 , I 
1250  SHORT  B(3),31,B2,33,B4,85,B6,F4,F5,F6,F9,F9,G5,37,3'3 
1260  SHORT  R0.R1 , R2 , R3 , R4 , R5 , R7 , R8 , T2 , T5 , T7 , T3 , U9 , V3 , V4 , 73 
1270  SHORT  X (50) , X 0 , Y ( 5 0 ) ,29 

1280  INTEGER  B8,89, 10, 13, 14, 15, 18, I9,K3,K9,N3,N4 , N9 

1290  INTEGER  SO, SI , S2 , S3 , S4 , S5 , S6 , S7 , S9 

1 3  0  0  1  initial  values 

1310  31 , S2 , S6= 1 

1320  33 , S4 , S9=0 

1330  D2~P!+i54#154/4  1  aaA2 

1340  G2 - 1  6  37 -28? . 639  9  G8  =  30 

1  350  R 0 s 5 1 2 . 5  9  Rl*600  0  R2*.6/R0  3  RS*1  3  T2=lS3  1  mo-e  c 

1  360  R  3  - . 0 1  !  F  a  d  j  u  s  t  error 

13“!'  R4, .95=1.01  '  peak  detect  error 

1330  14-6  9  GCSUB  2060 

1390  14-3  9  20=0  3  G0SU8  2060  1  position  input 
1400  14=2  9  GOGUB  2060  '  depth , temperature 
1410  R 7  = 1 20  '  3  outstep  ratio 

1420  B(i!,S(2!=2  9  B!3!=2.5  !  plot  scaling  parameters 
1430  15  =  5  9  T 9  =  0  1  read  delay  factors 
1440  K0=998.2  I  kg/m '3  water  density 


iaaaasaaaBsaaaasaassaaaB! 


!  begin:  equipment  setup33333333333333333333333333333— 

OUTPUT  709  |  "RS"  !  clear  A/D  unit 

OUTPUT  707  ; "PR"  !  clear  3314  fen  generator 

OUTPUT  723  j " H 2 "  !  xmit  voltmeter 

OUTPUT  707  ; "  AP10MV" 

OUTPUT  709  ; ” F1R0RA02 1N5LS5; T1 "  !  a/d  receiver  output 
OUTPUT  720  ; " F N 1 2 "  !  position  counter 
OUTPUT  709  ; "CL55" 

14  =  1  a  6QSUB  2060  !  LI  A  parameters 

!  end:  equipment  setup33333333333333333333333333333333' 

!  begin:  fresh  slate  for  processing  . . . 

13 , 19  ,  N9 , N4  =  0  !  plot  file  &wl/4  file  sizes  ,  wl / 4  count 

01,02,03,04,05,06=0  !  wavelength  sums 

F7  ,  F8 , 07 , 08 , 07 , V6=0 

Ml  ,112*1000000  I  0  dB  re  iV/uPa 

V 4  = .  0 0 0 1  9  V 3  =  1 0  9  T7  =  0 

S6= 1  !  ADC  on 

16=1  6  GOSUB  9120  !  switch 
GOSUB  1900  !  set  main  menu 

!  end:  fresh  slate  for  processing33-33333333--333"'--1 
!  begin  :  processing  stage3'333'3333333-'3-33333333-—333' 
!  hold  here  for  operator  menu  selection 
CLEAR 

D13P  USING  9590  ;  Al , F7 , F8 , 07 , 08 , 07 , FNM (Ml ) , FNM  ! M2 ' 

IF  S3  THEN  DISP  "plot  on  " ; ELSE  DISP  "plot  off 
IF  56  THEN  DISP  "ADC  on"  ELSE  DISP  "ADC  off" 

DISP  S*5m  Z0="iZ0 

KEY  LABEL 

FOR  17=1  TO  49999 

IF  S4  THEN  1760  !  fast  refresh  on  return  from  key  rout; 
NEXT  17 
S  4  -  0 

GOTO  1660 

!  end:  processing  stage3333--13  !-'-333 3333 3333-  »'‘ 3333  =  ’-  --  -  • 


fast  refresh  on  return  from  key  r  o  l 
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1800  !  begins  main  fflenua=3aaaaaaaaaaaaa= 

1810  ON  KEY#  1 ,  “  F  n "  SOSUB  2450 

1820  ON  KEY#  2,CHR*(11)  SOSUB  4150 

1830  ON  KEY#  3, "move"  GQSUB  7990 

1840  ON  KEY#  4,  "opr"  SOSUB  2060 

1850  ON  KEY#  5,"Q"  SOSUB  3380 

1860  ON  KEY#  6,"a/pt"  GOSUB  8120 

1870  ON  KEY#  7, "CAL"  GOSUB  5190 

1380  ON  KEY#  8, “clear"  GOTO  1930 

1890  S 4  =  1  !  fast  refresh' 

1900  RETURN 

1910  !  endsmain  menu'-aaaaaaaa  =  a  =  =  =  =  =  =  sa=< 

1920  !  begins  clear  slate  menu 

1930  !  provides  ex i t/ref resh  processing 

1940  OFF  TIMER#  1  0  OFF  TIMER#  2 

1950  OUTPUT  709  ;"QPNO"  S  OUTPUT  709  ;"0 

1960  S 1 , S2= 1  !  flag  reset 

1970  13,33,39*0  !  plot  damp 

1980  ON  KEY#  l,"main"  GOSUB  1800 

1990  ON  KEY#  5, "wipe"  GOTO  1550 

2000  ON  KEY#  6, "quit"  GOTO  8840 

2010  PRINT  "OPERATOR  CLEAR/RESET" 

2020  S 4 - 1  !  fast  refresh 

2030  GOTO  1650  !  enter  processing  stage 

2040  !  end:  clear  menu~“"=“ 


2400  !  title  (14=6) 

2410  A#*A*C1,183 

2420  DISP  "  d  a  t  e " ;  0  INPUT  W$  0  A*=A#&W*  0  GOSUB  8070 
2430  DISP  "run  #"  0  INPUT  W*@  A$=A«W*  0  GOSUB  8070 
2440  !  transducer  chars  only(i4=9) 

2450  At =A#C 1,643 

2460  H$  =  "T  is:"  ®  DISP  H* ; 0  INPUT  W$®  A#=AI&Hf&W*  S  GOSUB  807 
2470  DISP  "T  active  1  eng th ( mm ) , capaci tance (pF ) " ; ®  INPUT  D5,Y3 
2480  H$  = "  H  is:1'  0  DISP  H$;3  INPUT  W*0  A*  =  A*?<HJ &W* 

2490  DISP  "H  active  1 ength (mm) " ; Q  INPUT  D6 
2500  PRINT  USING  8630  ;  A*,D5,Y3,D6 

2510  Y3  =  Y3/10A12  0  D5=D5/ 1 000  0  D6  =  D6/ 1 000  I  SI  conversions 
2520  GOTO  2630 

2530  !  plot  enabl e/di  sab  1 e ( i 4  =  7) 

2540  S3=N0T  S3  0  13=0 

2530  IF  S3  THEN  2560  ELSE  2630 

2560  PLOTTER  IS  705 

2570  LIMIT  32,232,30,130  0  LOCATE  22,122,27,97  a  PEN  0 

2580  GOTO  2630 

2590  I  position(i4=8! 

2600  DISP  "position*"; ZO;0  INPUT  W$ 

2610  IF  LEN  (W$) >0  THEN  ZO-VAL(WI) 

2620  GOTO  2630  . 

2630  14=0  0  S4-1  0  RETURN 

2640  !  end:  operator  input  pr ocedur e=  =  =  =  =  =  =  =  =  =  =  =  =  =  = 5  =  =  =  ■  =  =  ■•  =  '; 


2650  !  begin:  peak  -frequency  measurement  pr ocedur 8=«a”"- 

2660  PRINT  USING  9460  j  A$  !  title 

2670  CLEAR  @  DISP  "peak  frequency  measurement" 

26B0  DISP  "esti mate , search  range, stepsize"; S  INPUT  F 7 , W 0 , F 9 
2690  F4  =  IP (F7-W0) 

2700  F5=CEIL (F7+W0) 

2710  !  documentation 

2720  PRINT  USING  9470  ;  T 1  ME , F7 , WO , F9 

2730  !  warmup  reading 

2740  IF  S5  THEN  V6  =  V3  ELSE  V6  =  Y4 

2750  F6=F7  @  T7  =  T8  0  GOSUB  6990  !  :<mtr  on  .long  delay 
2760  GOSUB  7310  !  warmup  read 

2770  IF  S3  THEN  S9=l  0  GOSUB  5770  !  begin  plot  of  a  vs  f 

2780  !  begin  sinitialize  sampling  array 

2790  F6=F4-3*F9  S  GOSUB  6990 

2900  FOR  12=1  TO  7 

2910  GOSUB  7310  !  read  amplitude  A2 

2820  IF  12(7  THEN  F6  =  F6  +  F9  3  GOSUB  6990  1  ;<mtr  on  next  fraq 
2830  A ( 1 2 ) =A2 
2940  NEXT  12 

2850  '  end  :fill  sample  array 
2860  !  begimresonance  search 

2870 -S7,F1 ,F2,F7  =  0  !  peak  detect  flag  and  detection  storage 
2330  FOR  F 0  =  F 4  TO  F5  STEP  F9 

2990  !  fO  is  the  freq  at  the  center  of  the  sample  data  set 
2900  F6=F0+4#F9  0  GOSUB  6990  !  xmtr  on  next  read  freq 
2910  !  curvature  at  fO 

2920  U2=(5*(A(7)+A(1> ) - 3  * ( A ( 5 ) +  A  ( 3 ) ) -  4  *  A ( 4 ) ) / 3  4 
2930  DISP  USING  8480  ;  F0,U2 
2940  IF  U 2 < 0  THEN  2950  ELSE  3120 
2950  1  slope  at  fO 

2960  Ul  =  (3*  (A  (7 ) -A  ( 1)  )  •►2*<A(6>  -  A  (2)  '  *  ( A  ( S  -A  (3)  j  !  /2S 
2970  !  peak  frequency 
2930  F3=F )-F9*U 1 / { 2 * U 2 ! 

2990  DISP  USING  8490  j  F3 

3000  IF  ADS (F3-F0) <F9*R4  THEN  3010  ELSE  3120 


!  peak  detected  at  center  of  sample  array 
!  peak  amplitude 

U0-(7*A(4)+6*(A(3>+A<5>)+3*(A(2)+A(6>)-2«<A(?>+A(U 
A7=U0-U1*U1/  < 4 * U 2 ) 

S7=S7+1  3  V7=Vt 

IF  S5  THEN  V3=V6  ELSE  V4»V6 

IF  S7= 1  THEN  F1=F3  ELSE  F2=F3 

PRINT  USING  8500  ;  F3, 10G0*A7*V1 t 1000*U2, 1000»Vl 
U9  =  F9  9  X0  =  F0 
GOTO  3160 

'  endipeak  detected  block 
!  begin:peak  not  detected  clock 
S7  =  Q  !  unset  peak  detect  flag 
1  end;  peak  not  detected 
!  end:tests  at  FO 
DISP 

IF  S7< 2  THEN  3180  ELSE  3240  !  end  search  on  two  det 
!  push  down  sampling  array 
FOR  10=1  TO  6 
A  ( 10) -A ( 10  +  1 ) 

NEXT  10 

G0SU8  7310  9  A(7)=A2 
NEXT  FO 
!  end  search 

IF  S3  THEN  G0SU3  5830  1  show  plot 

IF  S 7 < 2  THEN  COPY  <s  GOTO  3270  ELSE  3310 

IF  F1--0  THEN  PRINT  "search  failed"  9  GOTO  3350  ELSE 

!  one  good  reading 

F7-F1  3  F 3 3 F 9 

GOTO  3330 

1  two  good  readings 

F7=  (F1  +  F2)  /2  9  F9-»ABS  (F2-F1  >  '  average  of  two  good 

!  final  results 

PRINT  USING  9510  ;  F7.F9 

S 4 - 1  1  fast  refresh 

RETURN 

i  end:center  frag  mmt  procedure-s*—-  —  *3---3  -  =  = 


)  )  /  2 1 


!  begin!  Q  measurement  pr ocedurs==  =  =  *  =  ****!  =  =  =  s  =  =  =  =  =  “-  =  * 
CLEAR  3  DISP  "Q  measurement" 

DISP  "center  freq  ="|F7 
Q7=0 

PRINT  USING  S570  i  T IME , F7 , ZO , S* 

GOSUB  4950  !  measure  peak  voltage  v8  at  f7 
A7=*V0/Vi 

!  set  search  limits 
F9=CEIL (F7/R7) 

F4=MAX (0,F7-20*F9! 

F5*F7+204F9 

DISP  "peak  !  1000#A7*V1 

A4=sA7*.  707106791  185  !  1/2  power  amplitude 
DISP  "half  power  s " ; 1000*A4*V1 
!  begimsearch  for  lower  1/2  point 
F6 -  I P (F7) 

A8aA7 

!  begin  sloop  to  find  f  6<  f  1 
F6 -F6-F9 

IF  F 6 > F 4  THEN  GOSUB  6990  ELSE  4090  1  *.#tr  on,  or  e;<it 

GOSUS  7310  !  measure  A2 

DISP  "q  *;  ( A 2 - A 4 ) *1000*V1 

'  check  for  downhill  motion 

IF  A 2 > AS  THEN  3620  ELSE  3650 

1  uphill;  low  1/2  pt  obscured 

PRINT  USING  8520  ;  F6,A2,A8  6  COPY 

GOTO  4110 

!  downhi 1 1 -reset  check  val 
A8*A2 

IF  A20A4  THEN  3680  ELSE  3550 
'  end  loop 

'  begin  loop:  bracket  1/2  point 

IF  F6CF7  THEN  F6=F6+!  8  GOSUB  6990  ELSE  4090 

A8=A2  1  save  last  sample 

GOSUB  7310  !  new  sample  a2 

DISP  "q ! " ; ( A  2  -  A  4 ) *1000*01 

IF  A2 1-  =A4  THEN  3730  ELSE  3690 

1  end  loop. bracketed  -interpolate 

F1*F6-(A2— A4) / 

PRINT  USING  8530  ;  Fi 
'  end  search  for  FI 


."'V'  s'1.'. 
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3790  !  begin:  search  for  hi  1/2  power  pt  F2 
3300  F6=IP (F7) 

3310  A3*A7  I  check  value 
3820  !  begin:  loop  for  -f 6 > f 2 
3830  F6=F6+F9 

3340  IF  F 6 < F 5  THEN  GOSUB  6990  ELSE  4090 

3350  GOSUB  7310  !  measure  A2 

3360  DISP  "q  " ; ( A2-A4) *10Q0*V1 

3870  !  check  downhill  motion 

33S0  IF  A 2 > A 8  THEN  3990  ELSE  3920 

3390  !  uphill— hi  1/2  obscured 

3900  PRINT  USING  3520  ;  F6,A2,A8  3  COPY 

3910  GOTO  4110 

3920  !  downhi 1 1 -reset  check  val 
3930  A8  =  A2 

3940  IF  A20A4  THEN  3950  ELSE  3320 
3950  !  end  loop 

3960  !  begin  loop,  bracket  1/.2  point 

3970  IF  F6>F7  THEN  F6=F6-i  fi  GOSUB  6990  ELSE  4090 

3980  A8  =  A2  !  save  last  sample 

3990  GOSUB  7310  1  new  sample  a2 

4000  DISP  "q ! "; (A2-A4) *1000*01 

4010  IF  A2>»A4  THEN  4020  ELSE  3960 

4020  !  end  loop. brack  a ted-interpol ate 

4030  F2=F6+(A2*A4)/!A2-A8) 

4040  PRINT  USING  3540  :  F2 

4050  G7=F7/ (F2-F1 )  !  quality  factor' 

4060  S 4  =  i  !  fast  refresh 
4070  PRINT  USING  8550  ;  Q7 
4030  34=1  !  fast  refresh 
4090  RETURN 

4100  !  begin :error  processing 
4110  PRINT  USING  8560  ;  F6 
4120  COPY 
4130  GOTO  4080 

4140  1  end:Q  measurement  procedure**"***3”"*”"*5 
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4150  !  begin;  wavelength  measurement  pr ccedure====*==== 
4160  CLEAR  8  DISP  "min/max  positioner" 

4170  IF  S5  THEN  V6=V3  ELSE  V6  =  V4  !  select  T  or  P  drive 
4180  F6=F7  8  T7=TS  3  GOSUB  6990  !  xmtr  on, long  delay 
4190  GOSUB  7310  !  warmup  read 

4200  IF  S3  THEN  S9=2  3  GOSUB  5770  !  plot  start,  a  vs  z 

4210  DISP  "posi ti on=" ; Z0 

4220  DISP  "stepsi ze" ; 3  INPUT  Z9 

4230  PRINT  USING  8720  ;  T I  ME , Z0 , Z9 

4240  !  begin: initialize  sample  array 

4250  FOR  12=1  TO  7 

4260  GOSUB  7310  !  read  amplitude  A2 

4270  Z4*Z0  3  GOSUB  7610  I  move  down  z9  mm 

4280  Z (12) =  Z 4  3  A ( 12) =A2 

4290  NEXT  12 

4300  1  end:  initialize  array 
4310  S7=0 

4320  1  begin:  search  loop 
4330  Z4=Z(4!  !  array  center 
4340  !  parabolic  fit 

4350  U2=(5*iA(7)+A(l) )-3*(A(5)+A(3) )-4»A(4) )/84 
4360  Ul= (3* (A (7) -A(l) ) +2* (A (6) -A (2) )+(A(3)-A(3) ; ! /23 
4370  U0=(7*A(4)  *6*  (A  (5)  +  A  *  3 )  )  +3* (A (6! +A (2)  >-2*<Ai7)  ■>  ftil 
4330  U9  =  Z8  3  X 0  =  Z 4  3  A7  =  U0-U1*U1/ <4*U2) 

4390  Z2  =  Z4-Z3»U1/ ( 2 *U 2 )  !  predicted  optimum 

440)  DISP  USING  8660  ;  Z4,Z2,U1,U2 

4410  IF  87  OR  ABS ( Z2-Z4) <R3*Z9  THEN  4420  ELSE  4460 

4420  !  hit-  document 

4430  IF  S7  THEN  4550  1  end  search 

4440  S 7  = 1  1  set  search  end 

4450  Z3=Z2  !  save  peak  estimate 

4460  !  continue  search 

4470  '  push  down  sample  array 

443)  FOR  10=1  TO  6 

4490  Z ( 10) =Z ( 10+1 )  3  A ( 10) =A ( IC+l ) 

4500  NEXT  10 

4310  30SUB  7310  '  read  a2(z0) 

4520  2 ( 7 ) =  Z 0  3  A ( 7 ) =A2 

4530  IF  37  THEN  PRINT  USING  8730  ;  Z2 , 1 000*A7*V 1 , U 1 , U2 
7610 

4540  GOTO  4320 

4550  1  end:  search  loop 
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4560  !  print  final  estimate 

4570  PRINT  USING  8730  i  Z2 , 1000*A7*V 1 , U 1 , U2 

4530  22* (22+23) 12  !  average 

4590  PRINT  "average8" ; Z2j "mm" 

4600  IF  S3  THEN  GOSUB  5830  !  show  plot 

4610  !  begin  data  review 

4620  IF  S6  THEN  4630  ELSE  4640 

4630  IF  35  THEN  V38V6  EL.SE  V4=V6  1  store  ADC  results 
4640  DISP  CHR$ ( 1 1 ) ; " /4  count8" ; N9+2; S  INPUT  W$ 

4650  IF  LEN<W$)>0  THEN  4660  ELSE  N9=N9+2  0  GOTO  4680 
4660  W0=VAL (Wl) 

4670  IF  WOO  THEN  N9  =  W0  ELSE  4910 
4680  1  begimaccept  measurement 
4690  1 9 ~ 1 9  + 1 

4700  Mi ,L9-Z2/1000  !  meters 

4710  !  begins  sums 

4720  Cl -CON?  '  -x 

4 730  C2=C2  +  N9*N9  !  "'::A2 

4740  C38C3+L9  !  *y 

4750  C4=C4+N9*L9  1  '-tty 

4760  CS-«C5  +  L9*L9  1  -yA2 

4770  C6=I9*C2-Ci*Ci  !  del  factor 

4730  !  end;  sums 

4790  1  beg  i  h  :  c output e  sound  spd 

4200  IF  C 6 > 0  THEN  4010  ELSE  4370 

4510  1  two  points  c  r  more 

4320  W 1  - ( ! 9  *C4-C 1 *C3) / C  6  1  slope 

4530  NO" ( C2+C3-C 1 ♦ C A ) / C 6  1  intercept 

42  40  C 7 " 4 * M 1  *  F  7 

4530  IF  19/2  THEN  4860  ELSE  4870 

4S60  C3-4»50R(  (C5  +  1 9»W0  ♦  WO-*  H  l  *14  1  ♦C2  -2<N0*33  -2*W1  -C4*2  +  W0» 
19-2! #r?/C6»*F7 
4 2  7 0  1  document 

4  £  3 ;  PRINT  USING  369  0  ;  1 9 , C7 , C3 , 4000 *W1 
4590  CQTQ  4930 

“900  !  and  :  accept  measuremt 
■*9!0  1  begin:  discard  measuremt 
4920  PRINT  "mat  discarded" 

’930  S 4  - 1  0  RETURN 

4940  1  end:  wavelength  measurement  procedure-*8888-888 


4950  !  begin:  5-pt  peak  voltage  measurement  procedure 
4960  IF  S5  THEN  V6  =  V 3  ELSE  V6  =  V4  !  select  T  or  P  drive 
4970  F6=F7  9  T7=T8  9  GQSUB  6990  !  xmtr  on  ,long  delay 
4980  GOSUB  7310  !  warmup  read 

4990  !  beg i n : f i ve-poi n t  mmt  of  peak  amplitude  v2 

5000  S0  =  S6  9  S6  =  0  !  ADC  off 

5010  F6=F7-2  8  GOSUB  6990  !  set  xmtr  Fn-2 

5020  GOSUB  7310  !  read  v2 (Fn-2) ; 1st  term 

5030  F6=F6+ 1  9  GOSUB  6990  !  xmtr  on  Fn-t 

5040  VB—(3»V2/35>  '  first  term 

5050  GOSUB  7310  !  read  at  Fn-1 

5060  F6-F6+1  9  GOSUB  6990  !  xmtr  on  Fn 

5070  V8=V8+l2*V2/35  !  2d  term 

5080  GOSUB  7310  !  read  at  Fn 

5090  F6=F6+ 1  9  GQSUB  6990  !  xmtr  on  Fn  +  l 

5100  V8»V8+17»V2/35  1  center  term 

5110  GOSUB  7310  !  read  at  Fn+1 

3120  F  6  a  F  6  + 1  9  GOSUB  6990  !  xmtr  on  F n  +  2 

3130  V8»VB+l2*V2/35  !  4th  term 

5140  GOSUB  7310  !  read  at  Fn+2 

5150  V8-V8-3*V2/35  '  result  at  (7 

5160  36aS0  !  restore  ADC 

5170  RETURN 

5130  !  end:  peak  voltage  mmt . . . . . 


5190  !  begin:  calibrate  hydrophone  pr ocedur e========= 

5200  W 1 =C7/F7  !  wavelength 

5210  H 1  =  F N H (PI  *D5/W1 )  !  length  average  corrn  tor  T 
5220  H2=FNH  £  P I #D6/W  1 )  !  length  average  corrn  for  H 
5230  PRINT  USING  8700  ;  A*,FML<H11 ,FNL(H2) 

5240  !  begin: preamp  calibration 
5250  CLEAR  8  DISP  "calibration" 

5260  19=4  8  S0SUB  9120  !  switch 

5270  DISP  "preamp  calibrate"  8  WAIT  T8*3 

5230  F 6 3 F 7  8  V6  =  V5*38/32.  8  T7  =  TB  8  S05UB  6990  '  xatr 

5290  G0SUB  7310  !  read  L I A  out 

5300  G2=G7*A2 

5310  PRINT  USING  8740  ;  FNL (G2 ) , 1000*0  1 

5320  !  end:  preamp  calibration 

5330  !  TON  measurements 

5340  18=3  8  30SUB  8120  !  switch 

5350  DISP  " T > H "  8  INPUT  14$ 

5360  GQSUB  4950  '  5-pt  mint  of  v8 

5370  E3=V8*H2/G2  !  Eth 

5380  J3=V1*2*PI+F7#Y3  !  amperes 

5390  PRINT  USING  8750  ;  VI , 1000*J3, 1000000*E3 

5400  1  end  T1H 

5410  1  begin  P’H 

5420  1 8  =  -2  8  G0SIJ8  8120 

5430  DISP  " P 1 H "  0  INPUT  W$ 

5440  1  P  drive  adjust 

5450  F6  =  F?  9  V6--V4  8  T7=>T9  3  GOS'JB  6990  I  amtr  on 
5460  GQSUB  7310  1  warmup  reading 
5470  1  begin:  adjust  loop 
54  30  GQSUB  6990  !  :<mtr  on 

5490  GOS'JB  7310  8  E2*V2*H2/32  '  read  v2,  correct  for 
5500  DISP  USING  3710  ;  1000*V6, 1000*V2/62, l0000t0*E3 
3510  IF  A33IE2-E3> >R3#£3  THEN  V6  =  FNR  ! V6 »E3/E2 )  8  GOTO 
3320  '  end  adjust  loop 
3330  V 4 - V 6 

354  )  PRINT  "P  drive: " : V 4 ; " V " 

5 5 5 0  r  end  :  F  drive  adjust 

3560  30 SUB  4950  1  v8  mmt 

5570  E2=Ve»H2/G2  1  Eph 

5330  PRINT  "  Eph  =  "; 1000000*E2: "uV  ' 

5590  1  er:d  P1H 
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!  begin  P > T 

DISP  " P } T "  3  18  —  1  3  SOSUB  S 120  !  switch 

SOSUB  4950  !  vS  mmt 

E1=V9*H1/G2  !  Ept 

PRINT  "  Ept  ■  " j 1000000*E1 j "uV" 

!  end  PIT  mints 
!  begin  sensitivity  calc 
07=P  I  »D3*F7 / (K0*C7*C7*Q7) 

PRINT  USING  S760  ;  F7,C7,Q7 
PRINT  USING  9770  ;  D3,K0,J7 
M1»SBR(J7#E1*E3/(E2#J3) ) 
N2=SQR(J7*E2*£3/(E1*J3> > 

PRINT  USING  3780  ;  FNM  (Ml ) , FNM (H2) 

1  end  sensitivity  comp 
S 4  =  1  @  RETURN 

!  end:  cal  i  br  ate  phones  =  "  =  =  =="  =  — ==="  = 


5760  !  ===»ssssPLQT  PR0CSaaaaaaa 

5770  !  begin:  plot  initialize  pr ocedur eaaa=a=aaaa 

5780  DISP  "load  plotter"j9  INPUT  W*  1  dummy  input 

5790  IF  LEN ( W$ ) >0  THEN  S3a0  9  RETURN  !  abort  plot 

5800  PEN  1  9  FRAME  9  PENUP 

5810  RETURN 

5820  !  end:  plot  initialize  procedure33333"3*333 

5830  !  begin:  show  plot  procedure"333333"33"3 

5840  !  begin:  scaling 

3850  B l , B4=99999  9  82,85=-Bl 

5860  FOR  10=1  TO  13 

5870  W03X(I0)  9  W 1  =  Y ( 1 0 ) 

5830  IF  WO >82  THEN  B2  =  W0 
5890  IF  W0< B 1  THEN  E 1 =W0 
5900  IF  M 1 >85  THEN  35aWl 
5910  IF  W l < B4  THEN  B 4-W 1 
5920  NEXT  10 

5930  DISP  "data  range  x:";Bl;B2 

5940  DISP  "  y: " ; B 4 ; B 5 

5950  !  X  axis  scale 

5960  B3=5  !  x-interval 

5970  W 0 ” 5 * B 3  !  major  divisions 

5980  B8= (CEIL (B2/W0) -IP (B1/W0) )  *5 

5990  IF  B 8 > 5 0  THEN  6010  ELSE  6000 

6000  IF  B8< 20  THEN  6080  ELSE  6160 

6010  !  too  many  divisions 

6020  K3  =  1 

6030  B3=B3*B(K3)  0  K3*FNK(K3) 

6040  W0=5*83 

6050  B8=(CEIL(B2/W0)-IP(B1/W0) ) *5 
6060  IF  8 8 > 5 0  THEN  6030  ELSE  6160 
6070  !  end  too  many  divisions 
6080  1  begin  too  -few  divisions 
6090  K 3 •- 3 

6100  B3-B3/B(K3)  0  K  3  3  F  N  H  ( K  3 ! 

6110  IF  83-0  THEN  6120  ELSE  6200 
6120  >10=5*83 

6130  88= (CEIL ( B2/W0 ) - I P ( 3 1 /WO) !*5 
6140  IF  B8< 20  THEN  6100  ELSE  6160 
6150  !  end  too  tew  divisions 
6160  1  end  if 

6170  D i 3  IP (3 1  /WO)  *W0  !  low  :<  bdry 
6180  B2=B1+B3*B3  !  hi  :<  bdry 
6190  SOTO  6250 

6200  !  method  failed-set  default  x  scaling 
6210  B3-  ( B2-B 10  / 1 0 
6220  B9310 


DISP  "default  x-scaling" 

!  end  default  scaling 
!  end  x  scaling 
!  begin  y  scaling 
B6  = . 002  !  mV 
W2=5*B6 

B9  =  (CEIL (B5/W2) -IP (B4/W2) )  #5 
IF  B9>40  THEN  6320  ELSE  6310 
IF  B9C10  THEN  6380  ELSE  6450 
!  begin  too  many  divisions 
K3  =  3 

B 6 - B 6 * B (K3)  §  W2=5*B6  a  K3=FNK (K3) 
B9*(CEIL (B5/W2) -IP (B4/W2) )*5 
IF  B9 >40  THEN  6340  ELSE  6450 
!  end  too  many  divisions 
!  begin  too  few  divisions 
K3*i 

B6=B6/B  (K3)  a  W2=5*B6  3  K3  =  FNK  (K3) 
IF  B6 >0  THEN  6420  ELSE  6490 
B9= (CEIL  CB5/W2) -IP (B4/W2) >  *5 
IF  B?< 1 0  THEN  6400  ELSE  6450 
!  end  too  few  divisions 
!  end  if 

B4=IP (B4/W2) *W2 
B5=B4+86»B9  !  hi  y  bdry 
GOTO  6540 

'  method  fail-set  default  y-scalin 

86* ( B3-B4 ) /  1 0 

89-10 

DISP  "default  y-scaling" 

!  end  default  y  scaling 
!  end  y-scaling 
SCALE  B1,B2,B4,85 
1  end  scaling 
!  plot  points 
OUTPUT  705  ; "  S  M  o ;  " 

K3= 1 3-7 

IF  K3< 1  THEN  K3  =  0  9  GOTO  6640 
F0F  10=1  TO  K3 
PLOT  X  <  1 0 )  ,  Y ( 1 0 ) , 0 
NEXT  10 

OUTPUT  705  ;“SM*;'‘ 

K3= K3+1 

FOR  1 0  =  K3  TO  13 
PLOT  X  (10)  ,  Y  ( 1 0 )  ,0 
NEXT  10 

OUTPUT  705  j " S M ; " 

CSI2E  4  0  FXD  1,3 
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6710  LAXES  B3 , B6 , B 1 , B4 , 5 , 5 

6720  N 4  =  N 4  + 1  9  H*=A*C33,37]  9  PRINT  "plot  ";H*;"#";N4 
6730  SETGU 

6740  MOVE  0,62  9  LOIR  90  9  LQRG  6 

6750  LABEL  USING  ”K"  i  " amp  1 i tude (mV  rms>" 

6760  LOIR  0  9  LORG  4  9  HOVE  72,12 

6770  IF  S9*2  THEN  Wl="depth  (mm) "  ELSE  Wl ■ " f r equency ( Hz ) " 

6780  LABEL  USING  "K"  j  Ml 
6790  MOVE  72,5  9  CSIZE  5 

6800  IF  59-2  THEN  LABEL  USING  8670  ;  F6,SJ,HI,N4  ELSE  LABE-  U 
e680  |  ZO , SI , HI , N4 
6810  SETUU  9  PENUP 
6820  !  end:  point  plot 
6830  !  begin:  curve  fit  plot 
6840  B3»  <  B2-B1 ) / 100 
6850  FOR  W 1 =B 1  TO  B2  STEP ' B3 
6860  W2=W1-X0 

6870  W3a1000*Vl*(U0+Ul#W2/U9+U2*(W2*W2)/(U9*U9)) 

6880  IF  W3 >84  AND  W3<B5  THEN  PLOT  W1,W3 
6290  NEXT  Ml 
6900  PENUP 

6910  MOVE  0,0  9  PEN  0 
6920  !  end:  curve  Tit  plot 
6930  !  begin:  clear  plot 
6940  13*0  @  S9 -0 
6950  1  end:  clear  plot 
6960  RETURN 

6970  !  end:  plot  procedureaaaa==a=a===aa=aa==-aaaa----=aa 


69B0  !  =*====MEASUREMENT  PR0CEDURES=«””*””*=«“=*======»” 

6990  !  begins  amplitude  read  setup  pr ocedur  ea  =  s  =  *aa====s===="  - 

7000  !  input: 56,^6, V6  outs  timer  #1  set  for  read  enable,  read  en 
able  si  off 

7010  IF  S6  THEN  7020  ELSE  7250 
7020  !  begins  ADC  block 
7030  !  automatic  drive  adjust 

7040  OUTPUT  707  USING  8800  }  06 , F6  !  HP3314  drive  on 
7050  SI *0  8  ON  TIMER#  1,T7+I5*T5  GQSUB  7550  !  read  disable 
7060  !  begins  loop  to  adjust  V6  until  LIA  output  is  in  limits 
7070  WAIT  2 * T 5 

7080  ENTER  709  ;  A5  !  read  LIA  magni tude-i gnore  read  disable 
7090  IF  A5>. 8  THEN  7100  ELSE  7130 
7100  !  LIA  out  too  high 

7110  IF  ASM. 12  THEN  V6=FNR  (V6*.  5)  ELSE  VA-FNR  (V6*.  7) 

7120  GOTO  7170 

7130  IF  A5< . 1  THEN  7150  ELSE  7290  !  here's  the  loop  exst  point! 
7140  !  LIA  out  too  low 
7150  V6SMIN(10,FNR( 06*15) ) 

7160  GOTO  7170 

7170  1  reset  HP3314  and  retest 

7180  OUTPUT  707  USING  Se 1 0  ;  V6 

7190  OFF  TIMERS  1  3  Sl-0 

7200  ON  TIMER#  1 , T8+ ( 15-2! *T5  SOSUB  7550 

7210  PRINT  SIC1, 13)"  dr i ve" ; V6; "V" 

7220  !  hold  for  read  enable 
7230  IF  SI  THEN  7060  ELSE  7220 
7240  1  ends  ADC  block 
7250  !  begin:  no  ADC  block 
7260  OUTPUT  707  USING  8800  ;  V 6 , F 6 

7270  S i *0  3  ON  TIMER#  1,T7+I5»T5  GOSUB  7550  !  read  disable 
7200  !  endsno  ADC  block 
7290  T7 *0  0  RETURN 

7300  1  end:  amplitude  read  setup  procedure  »»«»»«!•»»»  » s 


7310  !  begin:  ampl  i  tude  read  pr  ocedur e“  =  =  =  =  s="  =  =  =  =  =  s==  =  =  = 

7320  !  ins  SI , VO  out:  v2,a2 

7330  !  hold  for  read  enable 

7340  IF  SI  THEN  7350  ELSE  7330 

7350  !  read  enable 

7360  TRIGGER  723  !  measure  xiatr  voltage 
7370  ENTER  709  }  A5  !  LIA  magnitude  output 
7390  V2»A5*V5  !  volts  receiver  output 
7390  ENTER  723  j  VI  !  xmtr  voltage 
7400  A2=V2/V1 

7410  !  check  for  LIA  overload 
7420  IF  A5>1  THEN  7430  ELSE  7450 
7430  !  alarm  -overload 

7440  PRINT  "LIA  OVERLOAD"  9  PRINT  F6;A5;V1  0  BEEP  55,2000 
7450  !  and:  overload  test 

7460  DI3P  USING  9590  ;  F6 , ZO , 1 000*V1 , 1000»V2 
7470  IF  S3  AND  S9  AND  I3< 50  THEN  7480  ELSE  7530 
7480  (graphics  storage 
7490  13=13+1 

7500  IF  S932  THEN  X  ( 13) =Z0  ELSE  X ( 1 3 ) =Fa 
7510  V ( 13! *1000*V2  !  mV 
7520  !  end  graphics  storage 
7530  RETURN 

7540  !  end  :  amplitude  read  procedure3333333333333333"" 
7550  !  bc-gin:read  enable  procedure3"3"""3333""” "-*■ 
7560  (  called  by  timerlfl-ressts  read  enable  flag  si 
7570  OFF  TIMER#  1 
7590  31=1 
7590  RETURN 

7600  !  end:  read  enable  procedure333333  —  333""333""3"3 


7610  !  begin!  move  manipulator  down  Z9  mm  =  3*3  =  32  =  333=="32  =  =  -"2 
7620  !  input!  z9 

7630  ENTER  720  j  NO  !  counter  re* 

7640  IF  ABS ( Z9 ) < . 5  THEN  76S0 

7650  IF  Z 9 > 0  THEN  OUTPUT  709  ;"CLSr'  §  N3»l  ELSE  OUTPUT  709  ;"CLS 
0"  3  N3*-l  3  Z9“-  ( Z9/R8) 

7660  ON  TIMER#  2 , IP ( Z9/R2-T2)  GQSUB  7700  !  set  stop  time 
7670  S1,S2*0  !  read  o*f,inmotion 
7690  RETURN 

7690  !  end:  move  manipulator  procedure333333333333333333"3333  13 
7700  1  begin:  manipulator  stop  procedure3222222233223"33"2"1-" 
7710  !  input  z0,z9,n0,n3 
7720  !  output  revised  zO,step  zB 

7730  IF  N3 >0  THEN  OUTPUT  709  ;“0PN1"  ELSE  OUTPUT  709  ;"0PN0" 

7740  OFF  TIMER#  2 

7750  ENTER  720  ;  N1 

7760  !  loop  until  ctr  stops 

7770  ENTER  720  ;  N2 

77S0  IF  N2=N1  THEN  7790  ELSE  N13N2  3  GDTO  7760 
7790  !  motion  stop 
7800  Z8=N3» (Nl-NO) / RO 
7810  Z0=Z0+Z8 

7920  DISP  USING  3790  ;  Z 0 , Z 8 
7830  S2 ~0  !  flag  motion  stop 
7340  1  auto  drive  control 
7350  IF  36  THEN  7860  ELSE  7930 
7360  !  begin:adjust  drive 
7870  ENTER  709  ;  A5 
7300  IF  A5 > . 95  THEN  7890  ELSE  7950 
7890  !  L I A  out  too  high 

7900  IF  ASM.  12  THEN  V6---FNR  ( V6 » .  5 )  ELSE  V6  =  FNR  (Vfi*.  7) 

7910  OUTPUT  707  USING  8810  ;  V6 
7920  PRINT  "dri ve: V6; "V" 

7930  WAIT  1 5*T5 *T3 
■’940  GOTO  7370 

7950  ON  TIMER#  1,T7+I5*T5  GOSUB  7530  !  set  delay  to  read  enable 
7960  77*0 
7970  RETURN 

7990  !  end : man i pu! at  or  stop  orocedure33"333333333"3333"113---13 
79-0  1  begin:  manipulator  relocate  pr  ocedure"""33""  -  ■  -  " " 

3  0  0  0  CLEAR 

3010  IF  Z0=0  THEN  DISP  "position"; 0  INPUT  ZO 

8020  DISP  "move  to"; 3  INPUT  Z4 

9030  Z9=Z4-Z0  1  mm 

8040  GQSUB  7610  !  move  :9  mm 

9050  S4= 1  3  RETURN 

8060  1  end:  manipulator  relocate  procedure33333333333'333"3"-13 
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8070  !  begin:  title  string  justifier====================== 

8090  1 0 = L EN (A$) 

8090  IF  10  MOD  32 >0  THEN  A$=A*&"  "  0  GOTO  8080 
8100  RETUPN 

8110  !  end  title  justif iar===============s====*====*====== 

8120  !  ADC/xmtr  toggle  routine==================-========= 

8130  IF  18=0  THEN  8140  ELSE  8170 

8140  !  operator  call 

8150  S 4  =  1  !  fast  refresh  menu 

8160  CLEAR  0  DISP  "1=PJT.  2=P1H  3=T1H  4=cal";0  INPUT  13 
8170  !  process  command 

8180  IF  1 9 > 0  THEN  S6»l  ELSE  S6=0  0  IB—IB  !  IS >0  turns  ADC 
8190  18=19  MOD  5 
8200  IF  18=0  THEN  9420 
8210  OUTPUT  707  ; "APOVQ" 

8220  ON  18  GOTO  8230,3230,8330,8330 
8230  !  PIT 

3240  S5=0  0  SJ="P>T" 

8250  OUTPUT  706  USING  8820  j  "B1A1" 

8260  OUTPUT  714  USING  8320  ;  ” 3 2 A 2 " 

8270  GOTO  8420 
8280  !  P>H 

8290  85=0  0  S I = " P } H " 

9300  OUTPUT  706  USING  8820  ;  "82A1" 

9310  OUTPUT  714  USING  8820  ;  "e2A2" 

9320  GOTO  8420 
8330  1  TJH 

3340  S5  =  1  0  SJ|="T>H" 

8350  OUTPUT  706  USING  3820  ;  "B2A2"  ! 

8360  OUTPUT  714  USING  8820  ;  " 9 3 A 3 ” 

8370  GOTO  8420 
3390  !  preamp  calibration 
8390  35  =  0  0  8 1  = " P > c a  1 " 

9400  OUTPUT  706  USING  8820  :  " B 3 A 1 " 

8410  OUTPUT  714  USING  8820  ;  " B7A 1 " 

9420  PRINT  USING  "  /  V. ! "  ;  Sf 
S430  13=0  0  RETURN 

94  40  !  and:  x m t r / a d c  toggle  routine  —  338  =  =  =  =  ===  =  =  =  =  =  =  = 


S450  !  Images  for  formatted  oatpat  =  =  !ia  =  s  =  =  =  =  ss  =  =  =  :s  =  =  =  =  ::  =  =  =!  =  "=:  =  “" 
8460  IMAGE  lllll K 

3470  IMAGE  / / " t i me: " , 60/ " cen t er  freq  measurement "/"estimate: ", 5D / 
"range: ",  30/ "step  si  os : ",30/ 

8430  IMAGE  40 , 2X , S . 4DE 

8490  IMAGE  #,80.20 

3500  IMAGE  /  "peak : " , 50. 20 , "Hz " /5X3D. 50 , "mV" /5XS. 3DE , " mV/Hz  A2" / " x m 
tr : " , K , "mV  “ 

3510  IMAGE  /"resonance  center  f requency"/5D. 20, ,20. 30, "H: " 
8520  IMAGE  /* 1/2  power  pt  obscured " /4D , 2 ( XK) 

8530  IMAGE  "low  1/2:", 50.0 

8540  IMAGE  "high  1/2:", 50.0 

3550  IMAGE  "resonance  quality  f actor :  " ,  20. 2D 

8560  IMAGE  /"search  limit  exceeded ",  50 

3570  IMAGE  /"time: " ,6D/"0  measurement"/"center  f req: " , 50. 20/"dopl 
h  " ,  4  0 . 0  /  K  / 

3580  IMAGE  " o" ,  40 ,  X4D. 0 , 60. 20 , 30. 50 

3590  IMAGE  IK!  "Fn  *" ,  K 20. 20/ "Cn-" ,  K ,"  */  " ,  20 .  20/ "  Dn  - l!  ,  K  ’ "  M  ' 
T '  a ,  K , "  d  B "  / "  M  ( H )  s  “ ,  K , "  d  B " 

3600  IMAGE  /"0-return  to  main  1 -L I  A" / "2'water  3” dr i v j " . " 

4-comments  3-constants" 

8610  IMAGE  "lockin  anal y: er " / " i nput  ", X2D. 40 , "mV" /" expand I 
:  me  ",  X5C  20  X50 ,  "ms " 

3620  IMAGE  "depth  (•«):“  ,K/"volu«e  (■'‘3!  0. 6D 

8630  IMAGE  K/"T  1  cr.g th  = " , K ,  "mm" / "capaci  tance -  " , K ,  "pF“ /  "H  length’-" 
■l,  "mm"/ 


2  b  4  0 

IMAGE 

"water  temper  at ure "/ 20 . 0 , "  C" 

2630 

IMAGE 

"6-titlo  7 -p l o t  cn/cff”/ 

"E-pcsi ti 

cn  5 “ 

dcr  specs" 

36a0 

IMAGE 

. , 4 0 . 2 0 , 2 X K / "  " , 2  ( XS . 3DE ) 

w  a 

IMAGE 

"  0  , "  H :  ",K,7X,K,"*",:D 

3  6  3':’ 

IMAGE 

30.0,  "  m  n  " ,  K  ,  6  X  ,  K , "  If "  ,20 

Q  *.  C  f*> 

IMAGE 

/  20,"  p  t  s "  /  "  sp  d ■ " , 4  0 . 20 , “  +  / - " 

'in 

20,  “ 

m/s"/ 

"  =  "  t-  "  ,r.  *  !l 

f  '  i  • 1 

3  ?  0  j 

IMAGE 

K / " T  carrn:",S2D.2D,"-JB"/"H  c 

o  r  r  n 

:",S 

n  n  *1  ^ 

1  -J  it  * 

87 1C 

IMAGE 

"pda  ”,  30. 20 , 2  C  X . 4DE ! 

3720 

IMAGE 

/ "  1 1  me : " ,  60  , "  ma::/min  locator 

H  /  H  ^ 

/  3 

tart 

:  !‘ ,  4D 

.  D,"  step:' 

:c.  0/ 

2  7  70 
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IMAGE 
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8820  IMAGE  #,4A 

8830  !  and:  1(83983==*“  ====  =  =  =  =  =  =  = 
8840  !  begin:  end  program  routine 
8850  OUTPUT  709  ; ” RS " 

8860  OUTPUT  707  ; "APOVO" 

8870  END 
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